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Abstract 



After introducing the perturbed metric tensor in a Friedmann— 
Robertson-Walker (FRW) background we review how the metric per- 
turbations can be spht into scalar, vector and tensor perturbations ac- 
cording to their transformation properties on spatial hypersurfaces and 
how these perturbations transform under small coordinate transforma- 
tions. This allows us to introduce the notion of gauge-invariant pertur- 
bations and to provide their governing equations. We then proceed to 
give the energy- and momentum-conservation equations and Einstein's 
field equations for a perturbed FRW metric for a single fluid and extend 
this set of equations to the multi-fluid case including energy transfer. 
In particular we are able to give an improved set of governing equations 
in terms of newly defined variables in the multi- fluid case. We also give 
the Klein-Gordon equations for multiple scalar fields. After introducing 
the notion of adiabatic and entropic perturbations we describe under 
what conditions curvature perturbations are conserved on large scales. 
We derive the "conservation law" for the curvature perturbation for the 
first time using only the energy conservation equation |105| ]. 

We then investigate the dynamics of assisted inflation. In this model 
an arbitrary number of scalar fields with exponential potentials evolve 
towards an inflationary scaling solution, even if each of the individual 
potentials is too steep to support inflation on its own. By choosing 
an appropriate rotation in field space we can write down for the first 
time explicitly the potential for the weighted mean field along the scal- 
ing solution and for fields orthogonal to it. This allows us to present 
analytic solutions describing homogeneous and inhomogeneous pertur- 
bations about the attractor solution without resorting to slow-roll ap- 
proximations. We discuss the curvature and isocurvature perturbation 
spectra produced from vacuum fluctuations during assisted inflation 

izi- 

Finally we investigate the recent claim that preheating after inflation 
may affect the amplitude of curvature perturbations on large scales, 
undermining the usual inflationary prediction. We analyze the simplest 
model of preheating analytically, and show that in linear perturbation 
theory the effect is negligible. The dominant effect is second-order in 
the field perturbation and we are able to show that this too is negligible, 
and hence conclude that preheating has no significant infiuence on large- 
scale perturbations in this model. We briefiy discuss the likelihood of 
an effect in other models ||6^ . 

We end this work with some concluding remarks, possible extensions 
and an outlook to future work. 
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1 Introduction 



It has almost become impossible to write about cosmology without having at 
least mentioned once that we are now experiencing and living through "the 
golden age" of cosmology. This thesis is no exception. 

The standard inflationary paradigm is an extremely successful model in ex- 
plaining observed structures in the Universe (see Refs. |6^, |6l|, ^ for reviews). 
In this cosmological model the universe is dominated at early times by a scalar 
field that gives rise to a period of accelerated expansion, which has become 
known as inflation. Inhomogeneities originate from the vacuum fluctuations 
of the infiaton field, which on being stretched to large scales become classical 
perturbations. The field inhomogeneities generate a perturbation in the curva- 
ture of uniform density hypersurfaces, and later on these inhomogeneities are 
inherited by matter and radiation when the infiaton field decays. Large scale 
structure then forms in the eras of radiation and matter domination through 
gravitational attraction of the seed fluctuations. 

Defect models of cosmic structure formation have for a long period been 
the only viable competitor for the inflation based models. In these models a 
network of defects, usually cosmic strings, seeds the growth of cosmic structure 
continuously through the history of the Universe. Unfortunately, at least for 
the researchers involved in the investigation of these models, recent cosmic 
microwave background (CMB) observations have nearly ruled these models 
out i. 

The amount of data available in cosmology is increasing rapidly. New 
galaxy surveys, like the 2-degree-Field and the Sloan-Digital-Sky-Survey are 
beginning to determine the matter power spectrum, the distribution of galaxies 
and clusters of galaxies, on large scales. New balloon based experiments, like 
Boomerang and Maxima, have already measured the fluctuations in the CMB 
with unprecedented accuracy. Forthcoming satellite missions, like MAP and 
Planck, are going to improve these results even more. 

Although it might seem that extending the parameter space of the cos- 
mological models that are investigated will lead to a decrease in accuracy, 
the quality and sheer quantity of the new data will nevertheless allow us to 
constrain the parameter space enough to rule out large classes of models and 
to give hints in which direction we should direct our theoretical and experi- 



mental attention in future |]13|. The combination of data from different CMB 
experiments with large scale structure and super-nova observations will be 
particularly useful in this respect. 

But to take full advantage from these new experiments and the new data 
they will provide one has to make accurate theoretical predictions. In this work 
we extend the theory of cosmological perturbations. We apply the methods we 
have developed to investigate some models of the early universe and calculate 
their observable consequences. 

The evolution and the dynamics of the universe are governed by Einstein's 
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theory of relativity, whose field equations are p2 



G^, = 87rGT^,, (1.1) 



where G^j, is the Einstein tensor, G Newton's constant and T^^ the energy- 
momentum tensor. From the contracted Bianchi identity G^,^!^ = 0, we find 
that the energy-momentum tensor is covariantly conserved. 



T,.: = 0- (1-2) 



In an unperturbed Friedmann- Robert son- Walker (FRW) universe, described 
by the line element 



ds"^ = -dt^ + a^{t) 



1 — nr 



:i.3) 



the homogeneous spatial hypersurfaces pick out a natural cosmic time coor- 
dinate, t, and hence a (3+1) decomposition of spacetime. Here a{t) is the 
scale factor and r, 9 and (f) are polar coordinates. But in the presence of 
inhomogeneities this choice of coordinates is no longer unambiguous. 

The need to clarify these ambiguities lead to the fiuid-fiow approach which 
uses the velocity field of the matter to define perturbed quantities orthogonal 



to the fiuid-fiow |3q, V% V3, 05 An alternative school has sought to 



define gauge-invariant perturbations in any coordinate system by constructing 
quantities that are explicitly invariant under first-order coordinate transforma- 



tions P, |5g, . Results obtained in the two formalisms can be difficult 

to compare. One stresses the virtue of invariance of metric perturbations un- 
der gauge transformations, while the other claims to be covariant and therefore 
manifestly gauge-invariant due to its physically transparent definition. 

We put special emphasis on the fact that perturbations defined on an un- 
ambiguous physical choice of hypersurface can always be written in a gauge- 
invariant manner. In the coordinate based formalism the choice of hy- 
persurface implies a particular choice of temporal gauge. By including the 
spatial gauge transformation from an arbitrary initial coordinate system all 
the metric perturbations can be given in an explicitly gauge-invariant form. 
In this language, linear perturbations in the fiuid-fiow or "covariant" approach 
appear as a particular gauge choice (the comoving orthogonal gauge 0, |55[| ) 
whose metric perturbations can be given in an explicitly gauge-invariant form. 
But there are many other possible choices of hypersurface, including the zero- 
shear (or longitudinal or conformal Newtonian) gauge |0, in which 
gauge-invariant quantities may be defined. 

Whereas there are still more things we do not know than we do know, and 
even more things we don't even know we don't know, it is safe to say that the 
age of "golden-ish" cosmology has begun. 
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2 Cosmological Perturbation Theory 



In this section we present the fundamental building blocks of perturbation 
theory in a cosmological context. After introducing the concept of gauge and 
gauge transformation we give the governing equations in gauge dependent and 
gauge invariant form for scalar, vector and tensor quantities in a perturbed 
Friedmann-Robertson- Walker (FRW) spacetime in the single fluid case. After 
a short section on scalar fields we then give the governing equations for multiple 
fluids. We conclude this section with an application of the theory to quantities 
that are conserved on large scales. 



2.1 Decomposing the metric tensor 

We consider first-order perturbations about a FRW background, so that the 
metric tensor can be split up as 

= gj^a + Sg^, , (2.4) 
where the background metric gj^J is given by 

C = oH„){-,\l), (2.5) 

and jij is the metric on the 3-dimensional space with constant curvature k. 
We will denote covariant spatial derivatives on this space as Xabc...\i Q- 

The metric tensor has 10 independent components in 4 dimensions. For lin- 
ear perturbations it turns out to be very useful to split the metric perturbation 
into different parts labelled scalar, vector or tensor according to their trans- 
formation properties on spatial hypersurfaces The reason for splitting 
the metric perturbation into these three types is that they are decoupled in 
the linear perturbation equations. 

Scalar perturbations can always be constructed from a scalar quantity, or its 
derivatives, and any background quantities such as the 3-metric jij. We can 
construct any first-order scalar metric perturbation in terms of four scalars (p, 
B, ip and E, where 

6goo = -a'20, (2.6) 
Sgoi = Sgio = a^B\i , (2.7) 
5gij = -2a\^ -iij - E\ij) . (2.8) 

Any 3- vector, such as B\i, constructed from a scalar is necessarily curl-free, 
i.e., = 0. Thus one can distinguish an intrinsically vector part of the met- 
ric perturbation Sg^i, which we denote by —Si, which gives a non- vanishing 
5go[i\j]. Similarly we can define a vector contribution to 5gij constructed from 

^See appendix for definitions and notation. 
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the (symmetric) derivative of a vector F{i\j). The scalar perturbations are dis- 
tinguished from the vector contribution by requiring that the vector part is 
divergence-free, i.e., 7*"'>S'j|j = 0. This means that there can be no contribution 
to Sqqo from vector perturbations to first-order. The decomposition of a vec- 
tor field into curl- and divergence-free parts in Euclidean space is known as 
Helmholtz's theorem. The curl-free and divergence-free parts are also called 
longitudinal and solenoidal, respectively. 

Finally there is a tensor contribution to 5gij = a?hij which is transverse 
{,l^^hij\k = 0), i.e., divergence-free, and trace-free {Y^hij = 0) which cannot 
be constructed from scalar or vector perturbations. 

We have introduced four scalar functions, two (spatial) vector valued func- 
tions with three components each, and a symmetric spatial tensor with six 
components. But these functions are subject to several constraints: hij is 
transverse and traceless, which contributes four constraints, Fi and Si are 
divergence-free, one constraint each. We are therefore finally left 10 new de- 
grees of freedom, the same number as the independent components of the 
perturbed metric. 

We can now write the most general metric perturbation to first-order as 



The contravariant metric tensor, including the perturbed part, follows by re- 

1-20) B,'-S' 



quiring (to first order), g^^, g^^ = 6/^, which gives 



9"' = a-'{v) 



fi| ' - (1 + 2ij)Y^ - 2E, ~F\' - F\ ' 




Thus we have the general linearly perturbed line element: 

+ [(1 - 2^)7i,- + 1E\ij + 2F,|, + dx'dx^ } . (2.11) 
2.2 Coordinate Transformation 



We use Eq. ( |2.11D as a starting point and perform a small coordinate trans- 



formation, where the new coordinate system is denoted by a tilde. 

The homogeneity of a FRW spacetime gives a natural choice of coordinates 
in the absence of perturbations. But in the presence of first-order perturba- 
tions we are free to make a first-order change in the coordinates, i.e., a gauge 
transformation, 

fi = n^i\ x^ = x^ + ei' + e, (2.12) 

where ^ = C,^{r],x^) and ^ = ^{ri,x'^) are arbitrary scalar functions, and ^* = 
^*(?7, X*) is a divergence-free 3- vector. The function determines the choice of 
constant-?] hypersurfaces, that is the time-slicing, while ^| * and ^* then select 
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the spatial coordinates within these hypersurfaces. The choice of coordinates is 
arbitrary to first-order and the definitions of the first-order metric and matter 
perturbations are thus gauge- dependent. 

The result of the transformation Eq. (|2.12|) acting on any quantity is that 
of taking the Lie derivative of the background value of that quantity, i.e. for 



every tensor T one has |101 



6T = ST + C^T^ 



; 



(2.13) 



where Tq is the background value of that quantity. 

We use a slightly different approach by directly perturbing the line element 
in Eq. ( p.ll| ). Starting with the total differentials of C,^ and ^ and c?^* 



we get using Eq. ( p.l2| ) 

drj = 
dx' = 



de = e'dfi+i\di' 

di = i'dff + i ijdx^ , 



df^-edf^-i\dx\ 

dx'-{^\' + C')dfi-{i\y+l\^)dx^ 



(2.14) 
(2.15) 
(2.16) 



where we have used the fact that ^°(?7,x*) = ,^°(r/,5;'') and C,{ri,x^] 

to first order in the coordinate transformations. 

Substituting the differentials Eq. (|2.17|) into Eq. (|2.11|) and using 



a{ri) = a{fj) - fa\ri) 



(2.17) 
(2.18) 



we then get for the line element in the "new" coordinate system, working only 
to first order in both the metric perturbations and the coordinate transforma- 
tions: 



ds' 



,2 / 



(77) {- (1 + 2 (0 - /^e" - e)) dv' + 2{B + e- dvdi' 

-2 {S, + el) df^dx^ + [(1 - 2 (V; + he)) 7^J + 2 (E - 0,,, 

+2 {Fiy - ^iij) + hij] di'dx^} . (2.19) 



The general form of the line element presented in Eq. ( 2.11 ) must be invariant 
under infinitesimal coordinate transformations, since the line element ds"^ is a 
scalar invariant. We can read off the transformation equations for the metric 
perturbations by writing down the line element in the "new" coordinate system 



ds' 



a iv) 



+ 



[1 + 24))dfi^ + 2{B\i - Si)dfjdx' 
[1 - 2^)7ij + 2Eiij + 2Fi\j + hij 



doc doc 



] ■ (2.20) 
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The coordinate transformation of Eq. (|2.12| ) induces a change in the scalar 
functions (j), ip, B and E defined by Eq. 



B 
E 



B+e-^', 



(2.21) 
(2.22) 
(2.23) 
(2.24) 



where h = a'/ a and a dash indicates differentiation with respect to conformal 
time 7]. The vector valued functions Si and Fi transform as 



Si 



(2.25) 
(2.26) 



The tensor part of the perturbation, hij, is unaffected by the gauge transfor- 
mations. 



2.3 Time slicing and spatial hyper surfaces 

We now relate the metric perturbations to physical, that is observable, quanti- 
ties [^. The unit time-like vector field orthogonal to constant-// hypersurfaces 
is 

N^ = -{l-(l),-B.' + S'), (2.27) 
and the covariant vector field is 

iV^ = -a(l + 0, 0) . (2.28) 

The covariant derivative of the time like unit vector field A^"^ can be decomposed 
uniquely as follows [|101|] : 

N^;u = a^u + ^9Pi.. - a^N, , (2.29) 
where the projection tensor P^j^y is given by 

P^. = g^. + N^N,. (2.30) 
The expansion rate is given as 

= N^.,,. (2.31) 

the shear is 

a,, = \p;^Pf {N^.,p + Np,^) - h)P,, , (2.32) 



Note that the covariant time-hke vector field defined in Eq. ( |2.2^ ) can not support 
vorticity. 
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and the acceleration is 

S = N^-.N'' . (2.33) 

Note that on spatial hypersurfaces the vorticity, the shear, the acceleration and 
the expansion coincide with their Newtonian counterparts in fluid dynamics 



83,111. 

The extrinsic curvature of hypersurfaces defined by A*"^ is defined as 

K^, = -P/N^-u ■ (2.34) 

The Raychauduri equation is given solely in terms of geometric quantities 
and gives the evolution of the expansion 6 with respect to proper time as 

— = a% - a^^a^^'- - -0' - R.^N^N'' . (2.35) 

UT O 

We can write the expansion, acceleration and shear of the vector field for a 
perturbed FRW spacetime, considering scalar quantities first, as 

e = 3^ (1 - 0) - -ip' - -V^ {B - E') , (2.36) 
ai = (2.37) 

aij = a (^a\ij - ^7ii , (2.38) 

where the scalar describing the shear is 

a = -B + E', (2.39) 

and where the time components are zero, i.e. Qq = 0, ctq^ = 0. Note that for 
an unperturbed background the expansion 6 coincides with the expansion rate 
of the spatial volume per unit proper time. 

The intrinsic spatial curvature on hypersurfaces of constant conformal time 
T] is given by [Q, ^ 

(3)^ = ^ + 1^^ + ±^2^ . (2.40) 

For a perturbation with comoving wavenumber k, such that V'^tp = —k'^ip, we 
therefore have 

5(3)/? = 4 (3/c - ^, (2.41) 

and ip is often simply referred to as the curvature perturbation. 

For vector perturbations we find = 0, aoo = o'oj = 0, and ^ = 0, and the 
only non-zero first order quantity is the shear given by 

r,, = a,, = { (5,,, + S',,,) + [F^j + F,,,)'} , (2.42) 

where we distinguished the vector quantities from the scalar ones by a "bar" . 

There is no tensor contribution to the expansion, the acceleration and the 
vorticity, but there is a non-zero contribution of the tensor pertubations to the 
shear, 

a., = ^a (2.43) 
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2.3.1 Evolution of the curvature perturbation 

We can now give a simple expression for the evolution of the curvature per- 
turbation ip. Multiplying Eq. ( ^.361 ) through by (1 + 0) in order to give the 



expansion with respect to coordinate time, t = / adrj, we get 

^ = (1 + 0)^ = 3if - 3^ + V^a , (2.44) 

where a = E — B /a. We can write this as an equation for the time evolution 
of ip in terms of the perturbed expansion, 66 = 6 — 3H, and the shear: 

^ = --S9 + -V^d. (2.45) 

3 3 

Note that this is independent of the field equations and follows simply from 
the geometry. It shows that on large scales (V^a 0) the change in the 
curvature perturbation, tp, is proportional to the change in the expansion 66. 

2.4 The stress-energy tensor for a fluid, 
including anisotropic stress 

Thus far we have concerned ourselves solely with the metric and its represen- 
tation under different choices of coordinates. However, in any non-vacuum 
spacetime we will also have matter fields to consider. Like the metric, the 
coordinate representation of these fields will also be gauge- dependent. 

The stress-energy tensor of a fiuid with density p, isotropic pressure p and 
4-velocity is given by (ESl 



n = {P + P) u^u, + p5^, + tt'^, , (2.46) 

where we have included the anisotropic stress tensor which decomposes into 
a trace- free scalar part, 11, a vector part, tt*, and a tensor part, (t™sor)^i^^ 
according to 

TT^ = n, ^ - + i (tT^i^. + TT^.f) + (t--^^);!^^. . (2.47) 

The anisotropic stress tensor has only spatial components, TTjj, and is gauge- 
invariant. The gauge- invariance can either be shown by direct calculation as 



in 1 55 1, or by observing that i^ij must be manifestly gauge-invariant due to 



its being zero in the background [^, since the background is FRW and so 
isotropic by definition. The linearly perturbed velocity can be written as 

= -[{l-(p), V.' + v"] , (2.48) 
Uy = a[-{l + (f)), v\^ + B\i + Vi- Si] , (2.49) 
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where we enforce the constraint ■u^m'^ = —1. We can introduce the velocity 

potential v since the flow is irrotational for scalar perturbations. We then get 
for the components of the stress energy tensor 

= -iPo + Sp), (2.50) 

T° = (po+Po) {B\, + v\, + v,-S,) , (2.51) 

n = -(Po+Po) (v^' + v^) , (2.52) 

= {p, + 6p)6% + n\. (2.53) 



Note, that the above definition differs slightly from that presented in Ref. 



Coordinate transformations affect the split between spatial and temporal 
components of the matter fields and so quantities like the density, pressure 
and 3- velocity are gauge-dependent, as described in Section |2.4.1| . Density 
and pressure are scalar quantities which transform as given in Eq. ( |2.54|) in 
the following section, but the velocity potential transforms a.s v = v + C,', as 
given in Eq. ( p.55| ). 



2.4.1 Transformations of scalar and vector matter quantities 

Any scalar p (including the fiuid density or pressure) which is homogeneous 
in the background FRW model can be written as p(?7, x*) = po{f]) + Sp{ri,x^). 
The perturbation in the scalar quantity then transforms as 

6p = 6p-ePo- (2.54) 

Physical scalars on the hypersurfaces, such as spatial curvature, acceleration, 
shear or the density perturbation Sp, only depend on the choice of temporal 
gauge, but are independent of the coordinates within the 3-dimensional 
hypersurfaces determined by ^. The spatial gauge, determined by ^, can only 
affect the components of 3-vectors or 3-tensors on the hypersurfaces but not 
3-scalars. 

Vector quantities that are derived from a potential, such as the velocity poten- 
tial V, only depend on the shift ^ within the hypersurface and are independent 
of We therefore find that the velocity potential transforms as 

v = v + ^'. (2.55) 

The function only affects the components of divergence-free 3-vectors and 
3-tensors within the 3-dimensional hypersurfaces, such as the velocity pertur- 
bation f*, which then transforms as 

= v' + C' ■ (2.56) 



2.5 Gauge-invariant combinations 

Scalar perturbations The gauge- dependence of the metric perturbations 
lead Bardeen to propose that only quantities that are explicitly gauge-invariant 
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under gauge transformations should be considered. The two scalar gauge func- 
tions allow two of the metric perturbations to be eliminated, implying that 
one should seek two remaining gauge-invariant combinations. By studying 
the transformation Eqs. (|2.21| - pl24D , Bardeen constructed two such quanti- 



ties g H I, 

$ = (j) + h{B - E') + {B - E'Y , (2.57) 
^ = i) -h{B - E') . (2.58) 

These turn out to coincide with the metric perturbations in a particular gauge, 
called variously the orthogonal zero-shear [0, ^ , conformal Newtonian |77 



or 



longitudinal gauge |8J]. It may therefore appear that this gauge is somehow 
preferred over other choices. However any unambiguous choice of time-slicing 
can be used to define explicitly gauge-invariant perturbations. The longitudi- 
nal gauge of Ref. [Bl] provides but one example, as we shall show in Section 



2.6. 



Vector perturbations There is one vector valued gauge transformation 
^* and we have introduced two vector functions Si and Fi into the metric, 
Eq. ( p^.lll ). From these we can construct only one variable independent of the 



vector valued functions which is given by 

S, = F, = S, + Fl. (2.59) 

2.6 Different time slicings 
2.6.1 Longitudinal gauge 

If we choose to work on spatial hypersurfaces with vanishing shear, we find 
from Eqs. ( p^.23D ,( p.24D and ( |2.39| ) that the shear scalar transforms as a = o"— ^'^ 



and this implies that starting from arbitrary coordinates we should perform a 
gauge-transformation 

= -B + E' . (2.60) 

This is sufficient to determine the scalar metric perturbations (j), ip, a or any 
other scalar quantity on these hypersurfaces. In addition, the longitudinal 
gauge is completely determined by the spatial gauge choice 

^i = E, (2.61) 

and hence E = B = 0. The remaining scalar metric perturbations 0, ip and 
the density perturbation 6p become 

0z = (f) + h{B-E') + {B-Ey, (2.62) 
= ^-hiB-E') , (2.63) 
6pi = 6p + p',{B-E') . (2.64) 



In Bardeen's notation these gauge-invariant perturbations are given as $ = $Ag(°) and 
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Note, that (pi and ipi are then identical to $ and \I/ defined in Eqs. ( p.57|) 
and (|2.58|) . These gauge-invariant quantities are simply a coordinate indepen- 
dent definition of the perturbations in the longitudinal gauge. Other specific 
gauge choices may equally he used to construct quantities that are manifestly 
gauge-invariant. 



2.6.2 Uniform curvature gauge 

An interesting alternative gauge choice, defined purely by local metric quan- 



tities is the uniform curvature gauge [^, 0, ^ |9^, also called the off- 
diagonal gauge |jl2|, |T^. In this gauge one selects spatial hypersurfaces on 
which the induced 3-metric is left unperturbed, which requires = E = Q. 
This corresponds to a gauge transformation 



h 



(2.65) 



The gauge-invariant definitions of the remaining metric degrees of freedom are 
then from Eqs. and (g^) 



(t) + ^ + 
B-E' - 



(2.66) 
(2.67) 



These gauge-invariant combinations were denoted A and B by Kodama and 
Sasaki 



. Perturbations of scalar quantities in this gauge become 

5pK 



(2.68) 



The shear perturbation in the uniform curvature gauge is just given by 
cr^ = —Bf^. This is closely related to the curvature perturbation in the zero- 
shear (longitudinal) gauge, ipi, given in equation ( p.63|) . 



B, = ~hi,i=C-il 



(2.69) 



Gauge-invariant quantities, such as or tfji are proportional to the displace- 
ment between two different choices of spatial hypersurface, which would vanish 
for a homogeneous model. 

In some circumstances it is actually more convenient to use the uniform- 
curvature gauge-invariant variables instead of $ and For instance, when 
calculating the evolution of perturbations during a collapsing "pre Big Bang" 
era the perturbations 0^ and 5^ may remain small even when $ and \1/ become 
large ^ 0. 

Note that the scalar field perturbation on uniform curvature hypersurfaces. 



(2.70) 



is the gauge-invariant scalar field perturbation used by Mukhanov 
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2.6.3 Synchronous gauge 



For comparison note that the synchronous gauge, defined by = 5 = 0, does 
not determine the time-shcing unambiguously. There is a residual gauge free- 
dom = X/a, where X{x^) is an arbitrary function of the spatial coordinates, 
and it is not possible to define gauge-invariant quantities in general using this 
gauge condition [^. This gauge was originally used by Lifshitz in his pio- 
neering work on perturbations in a FRW spacetime |6^. He dealt with the 
residual gauge freedom by eliminating the unphysical gauge modes through 
symmetry arguments. 



2.6.4 Comoving orthogonal gauge 



The comoving gauge is defined by choosing spatial coordinates such that the 
3-velocity of the fiuid vanishes, v = 0. Orthogonality of the constant-// hyper- 
surfaces to the 4-velocity, m^, then requires v + B = 0, which shows that the 
momentum vanishes as well. From Eqs. ( |2.23|) and (p.55|) this implies 



-{v + B), 

- / vdr] + i{x' 



(2.71) 



where ^(x*) represents a residual gauge freedom, corresponding to a constant 
shift of the spatial coordinates. All the 3-scalars like curvature, expansion, 
acceleration and shear are independent of C,- Applying the above transfor- 
mation from arbitrary coordinates, the scalar perturbations in the comoving 
orthogonal gauge can be written as 



<j) + ^[{v + B) a]' 
ip-h{v + B) , 
E+ vdr]~i. 



(2.72) 
(2.73) 
(2.74) 



Defined in this way, these combinations are gauge-invariant under transforma- 
tions of their component parts in exactly the same way as, for instance, $ and 
\l/ defined in Eqs. ( p.57|) and (|2.58|) , apart from the residual dependence of Em 
upon the choice of ^. 

Note that the curvature perturbation in the comoving gauge given above, 
Eq. ( p. 741) has been used first (with a constant pre-factor) by Lukash in 1980, 
|[70|| . It was later employed by Lyth and denoted 71 in his seminal paper 
and in many subsequent works, e.g. 0] and |64 



The density perturbation on the comoving orthogonal hypersurfaces is 
given by Eqs. ( |2.54| ) and (|2.71|) in gauge-invariant form as 



Spm = 5p + p'oiv + B) 



(2.75) 
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and corresponds to the gauge-invariant density perturbation emEoQ^^^ in the 
notation of Bardeen 0. The gauge- invariant scalar density perturbation A 



introduced in Refs. [|T^ and corresponds to ^p^^\l / Po- 

If we wish to write these quantities in terms of the metric perturbations 
rather than the velocity potential then we can use the Einstein equations, 
presented in Section ( p.8.3|) , to obtain 

h(f) + ip' - k(B - E') 

In particular we note that we can write the comoving curvature perturbation, 
given in Eq. (|2.74|) , in terms of the longitudinal gauge- invariant quantities as 

= P.77) 

which coincides (for k = 0) with the quantity denoted ( by Mukhanov, Feld- 
man and Brandenberger in 



2.6.5 Comoving total matter gauge 

The comoving total matter gauge extends the comoving orthogonal gauge from 
the single to the multi-fluid case. Whereas in the comoving gauge the fluid 
3-velocity and the momentum of the single fluid vanished, in the total matter 
gauge the total momentum vanishes, 

(p + p)(v + b)= J2iPa + Pa) {ya + B^=0, (2.78) 

a 

where Va, Pa and pa are the velocity, the density and the pressure of the fluid 
species a, respectively. Orthogonality of the constant-?] hypersurfaces to the 
total 4-momentum, u'^, then again requires that -8 = 0, independently. Note 
that the gauge-invariant velocity V = v + E' (actually the velocity in the 
longitudinal gauge) coincides with the shear a of the constant-?] hypersurfaces 
in the total matter gauge . 

2.6.6 Uniform density gauge 

Alternatively, we can use the matter content to pick out uniform density hy- 
persurfaces on which to define perturbed quantities. Using Eq. (|2.54| ) we see 
that this implies a gauge transformation 

a = ^- (2.79) 

Po 

On these hypersurfaces the gauge-invariant curvature perturbation is p2|, ^ 



-( = ^Sp = iJ + h^. (2.80) 
Po 

The sign is chosen to coincide with ( defined in Refs. ^. There is another 
degree of freedom inside the spatial hypersurfaces and we can choose either B, 
or {; to be zero. 
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2.7 Gauge-dependent equations of motion 



We now give the governing equations in the homogeneous background and 
the gauge-dependent equations of motions for the scalar, vector and tensor 
perturbations. In the following we use the stress energy tensor of a fluid, 
including anisotropic stresses, as given in Section [2^. 



2.7.1 Conservation of the energy- momentum tensor 

In the homogeneous background the energy conservation equation is given by 

p'=-3h{p + p). (2.81) 

Note, there is no zeroth order momentum conservation equation as momentum 
is zero by assumption of isotropy. 

The conservation of energy-momentum yields one evolution equation for 
the perturbation in the energy density 

6p' + 3h{6p + 6p) = {p + p) [3iP' -V^{v + E')] , (2.82) 

plus an evolution equation for the momentum 

[(p + p){v + B)\ + 5p+\{V^ + 3K)Yi = -{p + p)[<j) + Ah{v + B)] . (2.83) 

There is no equivalent to the energy conservation equation in the case of vector 
perturbations since energy is a scalar quantity, but Eq. ( |1.2| ) gives a momentum 
conservation equation for the vector perturbations, 

[{^p + p){v,-Si)'{ + Ah{p + p){v,-Si) = -V,(4 + 7r,|^) 

= -(V' + 2fi:)7r,, (2.84) 

since momentum is a vector. 

There is neither an energy nor a momentum conservation equation for the 
tensor matter variables, since, as already mentioned, energy is a scalar and 
momentum a vector quantity, and hence they are decoupled from the tensor 
quantities. 



2.7.2 Einstein's field equations 

The equations of motion for the homogeneous background with scale factor 
a{r]) and Hubble rate h/a = a'/a^ are 

= ^pa^-^^ (2.85) 

h' = _l!L^(p + 3p)a2_ (^2.86) 
3 
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They are derived from the — and the i — j components of the unperturbed 
Einstein equations, respectively. Equation ( p.85|) is often referred to as the 
"Friedmann equation" . 

We now give the perturbed Einstein field equations beginning with the 
scalar perturbations in the single fluid case. The first-order perturbed Einstein 
equations yield two evolution equations from the i—j component and its trace, 
respectively, 

ip" + 2hil)' -Kip + h(f)' + {2h' + h^)(p = 4nGa^ (dp + -V^H ) , (2.87) 



3 

a' + 2ha~^ + ip = SiiGa'^U, (2.88) 

plus the energy and momentum constraints 

3h{ip' + h(f))-{V^ + 3K)iP-hV'^a = -AirGa^dp, (2.89) 

ip' + h^ + na = -4:7cGa'^{p + p){v + B), (2.90) 

derived from the — and the — i components of the Einstein equations, 
respectively, where a = —B + as given in Eq. (p.39|). 



The — i component of the Einstein equation for vector perturbations gives 
rise to a single constraint equation, 

(V' + 2k) {Si + F'^) = 16nGa\p + p) {S, - v^) . (2.91) 

The i—j component then leads to an evolution equation, expressed in terms 
of the vector-shear Tij defined in Eq. ( |2.42| ), 

4 + hTij = a^2TTG (vr,,, + 7r,|i) . (2.92) 

The case where the matter content of the universe is dominated by a scalar 
field is of great interest in cosmology. In anticipation of the results of Section 
|2.9| we can now show that in this case the vector perturbations are zero and 
stay so. From Eqs. ( ^.51| ),( P3BD and ( ^.118| ) and ( ^.119| ) we see that the vector 
matter variables Vi — Si and VTj are zero and there are therefore no source terms 
in the constraint equation ( |2.91| ) and the evolution equation ( |2.92| ). Hence the 
vector perturbations are zero in all of space and stay so as long as the scalar 
field dominates the energy content of the universe and afterwards if no sources 
of vorticity are introduced. 

The only non-zero component of the perturbed Einstein tensor for tensor 
perturbations is SG^ j. This gives rise to 

[h'lj + 2h h[j + {2k - V^) hij] = IQnGa^^'^'^'^'^TCij . (2.93) 

There is no separate conservation equation for the tensor matter variables. 
Hence the time evolution of the tensor metric perturbations hij is determined 
solely by the Hubble parameter h and the scale factor a, only sourced by 
matter perturbations in form of the tensor-anisotropic stress ^^'^^^"'^^ Hij . For 
scalar and vector matter fields ^^'^^^°^^7Cij = for linear perturbations and the 
evolution equation for hij is source-free or homogeneous. 



15 



2.8 Picking a gauge: three examples 



To illustrate the gauge invariant approach we use the gauge dependent equa- 
tions derived in the last section and present them in terms of gauge-invariant 
quantities which coincide with physical quantities in particular time-slicings. 
As examples we work with quantities coinciding with three specific gauges: 
the popular longitudinal gauge, the uniform density gauge and the comoving 
gauge. 



2.8.1 Governing equations in the longitudinal gauge 



The longitudinal gauge is defined by vanishing shift vector, B = 0, and vanish- 
ing anisotropic potential , E = 0, and hence ai = 0. Note that the influential 
report by Mukhanov, Feldman and Brandenberger, employs this gauge 
throughout. We now give the gauge invariant equations of motion in this 
particular gauge. 

We get from the conservation of the energy momentum tensor a "continu- 
ity" and a constraint equation. 



Sp[ + 3h{5pi + 5pi) 
[(p + p)^i]' + 5pi + ^(V' + 3K)fli 
The evolution equations become 

< + 2htp[ - + h4>[ + {2h' + = 

V^i - 01 = 

and the energy and momentum constraints are 



-{p + p) 



+ 4:hv\ 



3h{^[ + h^i) - (V^ + ?>k)^i 



—AixGa^dpi , 
—ATxGa^{p + p)v\ . 



(2.94) 
(2.95) 

,(2.96) 
(2.97) 



(2.98) 
(2.99) 



From Eq. ( p.97| ) we see that in the case of vanishing anisotropic stresses, as 
is the case for a perfect fluid or a scalar field, the curvature perturbation and 
the lapse function coincide, 

4^i = 4>i. (2.100) 

This can simplify calculations considerably. 



2.8.2 Governing equations in the uniform density gauge 

In the uniform density gauge the density perturbation vanishes, 6p = 0. This 
defines the spatial hypersurfaces. There is another degree of freedom inside 
the spatial hypersurfaces and we can choose either B, E ot v being zero. We 
choose the shift vector B = 0, so that asp = Esp. We now give the gauge 
invariant equations of motion in this gauge. 
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The conservation of the stress energy tensor gives then two constraint equa- 
tions 



3h 
p + p 

[{P + P)v5p]' + SpSp 



3V'L-V {vsp + asp 



Sp 



-{p + p) 



5<5p 



+ 4:hv. 



5p 



(2.101) 



(2.102) 



It is worth to point out the importance of the perturbed energy conservation 
equation in the uniform density gauge, Eq. ( |2.101| ) above, for later parts of 
this work. In Section |2.12| we discuss conserved quantities on large scales. 
Postponing a detailed discussion to the later section, we can already see that 
on large scales, i.e. when V'^{vsp + a^p) —>■ 0, the change in the curvature per- 
turbation is proportional to the pressure perturbation, ipg^ oc Spsp. Hence the 

curvature perturbation ipsp is constant on large scales for a vanishing pressure 
perturbation. In Section ^ on preheating we will discuss the case in which the 
pressure perturbation does not vanish and might give rise to a change in the 
curvature perturbation on large scales. 

The Einstein evolution equations in the uniform density gauge are 



A-KGa' [ 6psp + - V'n^p 



a'sp + 2hasp - (psp + i'Sp = STrGa^II, 



sp ' 



and the constraint equations are 



3h{^'sp + h<f)sp)-{V^ + 3K)^Psp - hV^asp = 

ip'sp + h(j)sp + Kcrsp = -ATxGa^{p + p)v5p . 



(2.103) 
(2.104) 

(2.105) 
(2.106) 



2.8.3 Governing equations in the comoving gauge 

The comoving gauge is defined by vanishing fiuid 3-velocity, 5 = 0, and van- 
ishing scalar shift function, 5 = 0. 

We now give the gauge invariant equations of motion in this particular 
gauge. We get from the conservation of the energy momentum tensor a "con- 
tinuity" and a constraint equation, 



Vm + 3/l(5pm + 5Pra 

SPm + ^(V^ + 3k)U, 



(p + p) ^fm 
-{P+P)4>m- 



(2.107) 
(2.108) 



The evolution equations become 
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and the energy and momentum constraints are 



0, 



(2.109) 
(2.110) 



(2.111) 
(2.112) 



We will use an extension of the comoving orthogonal gauge, the total matter 
defined in Section 2.6.5 , in Section 2.10| on multi-component fiuids and will 
give the governing equations in the total matter gauge there. 



2.9 Scalar fields 

In this section we briefiy introduce the stress energy tensor for scalar fields, 
postponing a detailed discussion of scalar fields and their dynamics to Sec- 
tion Ij. 

2.9.1 Single scalar field 

A minimally coupled scalar field is specified by the Lagrangian density 

C = ~^-'^y,.,,-Viv), (2.113) 

where ip'^ = g^^ip^^. The energy momentum tensor is defined as 

dC 

T% = -2-—^b^L, (2.114) 
and we therefore get for a scalar field 

n = - 5^ {v{^) + ^y^^V;.) • (2.115) 

Splitting the scalar field into a homogeneous background field and a pertur- 
bation, 

^{'n,x') = ^o{v) + 5^{'n,x'), (2.116) 

and using the definitions above we find for the components of the energy 
momentum tensor of a perturbed scalar field without specifying a gauge yet 

n = -^a-Vo'-^o + a-VU^V'o- V)-V;^5¥', (2.117) 
T\ = -a-^^',5v,), (2.118) 



S',, (2.119) 
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where = and Vq = ^^(v^o) 0- We also see, by comparing Eq. ( p.ll9|) 



with Eq. ( |2.53| ) above, that scalar fields neither support anisotropic stresses 



nor source vector and tensor perturbations, to first order. 

The Klein-Gordon equation or scalar field equation can either be derived 
by using the conservation of the stress energy tensor, T'^'^^^ = 0, or by di- 
rectly varying the action for the scalar field, equation ( p.ll3| ). We find for the 
background field 

^'^ + 2hip'Q + aW^ = , (2.120) 
and for the perturbed Klein-Gordon equation for the field fiuctuation 

Sip" +2h5ip' -\/'^5ip+a^V^^5ip+2a'^V^(l)-iPo(l)' -ipo {Sf - VV) = . (2.121) 
2.9.2 scalar fields 

For minimally coupled scalar fields the Lagrangian density is given by 

'^ = -^E«^^;/^)-^(^^)- (2.122) 

For the general potential V = V{ipi) the energy momentum tensor can not be 
split into a background and a perturbed part. In the special case of an additive 
potential of the form V = ^^i^ possible, but we postpone this 

case until we investigate the assisted infiation model in Section ||. 
The total energy-momentum tensor in the background is hence 



n = -a~' 



= 0, (2.124) 
T] = a- (^5^ - avj 5V (2.125) 

The perturbed energy-momentum tensor can be be given for each field / and 
is 

^(7)0 = -a~'{^'joSy^'j-^%^ + aW^^6^j) , (2.126) 
T°)^. = -a-Vfo^^/,,, (2.127) 
Til) J = a-'tK^. W + V'/o (M - V'^o)] , (2.128) 

where V^j = As in the single field case there are no vector and tensor 
parts and no anisotropic stresses to first order. The total perturbed energy 



^Note that there is a typographical error in Eq. (6.6) of Ref. ||8J], in which Mukhanov, 
Feldman and Brandenberger define their scalar field energy momentum tensor. 



19 



momentum tensor for all the fields is then given as the sum of the perturbed 
energy momentum tensors of each field, 

6n = Y,STi;^.- (2.129) 
I 

The background Klein-Gordon equation for the J-th scalar field is 

^% + 2h^'jo + a V^, = . (2.130) 
The perturbed Klein-Gordon equation for the J-th scalar field is 

S<f'l +2h 6if'j - V^6ifi + a^Yl ^mvMj 

J 

+2aV^> - ^',00' - (^',0 (3^' - VV) = , (2.131) 
where V^^pj^pj = q^^^q^j and the equation is gauge dependent. 

2.9.3 The Klein-Gordon equation in specific gauges 



As in Section |2.8| , we give here the perturbed Klein-Gordon equation ( p.l21| ) in 
gauge-invariant form in three different gauges, the longitudinal, the comoving 
and the constant energy gauge. We give the equations for the case of a single 
scalar field, which simplifies the notation, but they can be readily extended 
to the multi-field case, as can be seen by comparing Sections 2.9.1 and 2.9.2 
above. 

The Klein-Gordon equation in the longitudinal gauge is 



+ 2h5y^ - V^5^^ + a^V^/^^ + 2a^V^^x - - ^^'^i)'^ = . (2. 132) 

The Klein- Gordon equation in the uniform density gauge can be simplified 
by using the fact that 5p and 5p vanish. We find therefore in addition to the 
Klein-Gordon equation two constraints, 

6~^l - 2h6^,p + a^V^/^,p-2^l<Psp-Mp = ^- (2-133) 

Hence the perturbed Klein-Gordon equation in the constant energy gauge re- 
duces, subject to Eqs. ( p.l33| ), to a constraint equation in the field fluctuation, 

Aha^VJ^Sp + v'^VH^sp - 2(^'oVo05p - {^i^'sp - W^p) = , (2.134) 

which is not very useful, really. 

In the comoving gauge we find by comparing Eqs. ( p.51| ) and ( 2.118 ) that 
the field fluctuations vanish, 

5^^ = Q. (2.135) 

The perturbed Klein-Gordon equation then reduces to a constraint equation 
for the metric perturbations, 

2a V^0^ - ^'^^'^ - i^ij'^ - V V^) = , (2.136) 

where a^ = E'^. 
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2.10 Multicomponent fluids 

2.10.1 Including energy and momentum transfer 

In this section we will consider the evolution of linear scalar perturbations 
about a homogeneous and isotropic (FRW) universe, now containing several 
interacting fluids. 

The metric perturbations are the same as in the single fluid case, whereas 
the matter variables are different in the multi-fluid case for each fluid. The 
Einstein equations can be taken from Section (|2.7| ) (for each fluid species), 
but the energy and momentum conservation equation acquire new interaction 
terms. 

The total energy momentum tensor is given as 

T^u = . (2.137) 

a 

As in the single fluid case the total energy momentum tensor is conserved, 

T'^';^ = 0, (2.138) 
but the (no n-) conservation of energy- momentum for each fluid leads to 

Tf:y,, = QU^ (2.139) 

where Q^^^) energy momentum four vector. We can now give the govern- 

ing equations in the background. The energy momentum four vector in the 
background is given as 

^°^^?W = (-«QW' 0.0,0). (2.140) 

The energy conservation equation for fluid a with energy density pa and pres- 
sure Pa is then 

p'a = -^h{l - qa)ipa + Pa) , (2.141) 

where 

parameterises the energy transfer between fluids, subject to the constraint 

5^(Pa+Pa)ga =0, (2.143) 

a 

which ensures that the total energy is conserved, as in the single fluid case. It 
follows from Eq. ( p.l37| ) that the total density and the total pressure are the 
sum of the densities and pressures, respectively, of each fluid component a, 

P = Y^Pa. P = Y^Pa. (2.144) 
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The equations of motion for the homogeneous background are the same as in 
the single component case, see Section (|2.7| ). The perturbed energy momentum 
four vector is given as 



(a) 



6Q- 



-aQ{a){(t) + Ca) , 



(2.145) 



where we introduced ea and /(a) to describe the energy and momentum trans- 
fer, respectively. Under a scalar gauge transformation transforms as 



/ 



a) 



(2.146) 



whereas the momentum transfer parameter /(„) is gauge invariant. 



2.10.2 Scalar perturbations in the multi-fluid case 

We get for the gauge dependent perturbed energy (non)-conservation equation 

6p'^ + 3h {6p^ + 6po,) = {pa + Pa) [- V (v^ + E') 

+3'?/'' + 3hqa{(j) + Cq 



(2.147) 



and for the perturbed momentum equation 

[{pa + Pa) (Va + B)]' + (p« + p«) [Ah{Va + B) + (f)] 

+5Pa + ^(V' + 3K) n„ = h{p^ + p^) [3q^{v + B) + F^] 



(2.148) 



Since the metric is the same for all the fluids the field equations do not change 
compared to the single fluid pointed out in the previous section. 

2.10.3 Vector perturbations in the multi-fluid case 

We extend the equation of motion for vector perturbations in the single fluid 
case, Eq. (|2.84|), to the multi- fluid case: 



[{pa+Pa)Vai]' + 4:h{pa+Pa)Vai-- (V^ + 2k) TC^i = Kp^+p^) fat + QaVa 

(2.149) 

where Vai = V(^a)i — Si is the gauge-invariant vector velocity perturbation, 
fai is the gauge-invariant vector momentum transfer, and iiai = 7i(^a)i is the 
anisotropic stress of the species a, respectively. 
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2.10.4 Tensor perturbations in the multi-fluid case 

The tensor perturbations are coupled directly to the matter only through the 
anisotropic tensor stress. In order to extend the equation of motion for tensor 
perturbations, Eq. (|2.93|) , to the multi-fluid case we therefore simply have to 



set 

(tensor) 



^u = E^*"'"^^W*^' (2.150) 



where ^^"^^^""^^7^(0)^ is the anisotropic stress perturbation for the tensor mode of 
species a. 

2.10.5 Equations of motion for scalar perturbations in the total 
matter gauge 

Although we work in the total matter gauge throughout this section we omit 
the "tilde" and the subscript "m" to denote the chosen gauge, since there is 
no confusion possible and the equations appear less cluttered. 

The first-order perturbed Einstein equations in the total matter gauge yield 
two evolution equations 

^" + 2hi)' + h(j)' + {2h' + h^)(t) = AtiGo^ (^p + ^V^n^ ,(2.151) 

a' + 2ha - (j) + = SnGa'^Il, (2.152) 

plus the energy and momentum constraints 

3h{ip' + h<p) - (V^ + 3k)V^ + WV = -A^Ga^dp , (2.153) 

tlj' + h(i)+Ka = 0. (2.154) 

The (non-) conservation of energy- momentum for each fluid leads to 

T(^a) 0-^, = -3/i(Pa+Pa)ga(l+0 + ea), (2.155) 
^(a) = h{Pa+Pa)fa\i, (2.156) 

SO that the perturbed energy transfer in the total matter gauge is given 3h{pa + 
PajQa^a and the momentum transfer is given by the gradient of h{pa +Pa)fa- 
This yields the equations of motion for the density and velocity perturba- 
tions 

6p^ + 3h {6pa + Spa) = {Pa + Pa) ["V^ {v^ + O") 

+ 3{tlj' + h(f)) + 3hq^ea - 3h{l - q^)(f)] ,(2.157) 

2 

[ipa+Pa)Va]' = (pa+Pa) [-Ahv^ - (p + hf^] - Spa - -{V^ + 3k)U^ . (2.158) 

The conservation of the total energy-momentum yields one evolution equation 
for the density perturbation 

6p' + 3h{p + p){6p + 6p) = (p + p)(3^'- VV), (2.159) 
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plus a momentum conservation equation, which in this gauge reduces to an 
equation for the hydrostatic equihbrium between the pressure gradient and 
the gravitational potential gradient, 



(2.160) 



Use of the energy and momentum constraint equations ( 2.153| ), ( p. 154 ) and 
( p.l60| ), considerably simplifies the evolution equations ( |2.159| ) and ( p.l52| ) to 



give 



6p' + 3h6p 
a' + ho 



-(V^ + 3k) [{p + p)a - 2hl{] 



+ mGa — 



p + p 



3 p + p 



(2.161) 
(2.162) 



where the gauge-invariant metric perturbation \1/ = ip + ha (actually the curva- 
ture perturbation in the longitudinal gauge, Eq. (|2.63|) is related to the density 
perturbation in the total matter gauge by the energy constraint equation ( p. 891) 



47rGa2 
F -3k 



5p. 



(2.163) 



2.10.6 Entropy perturbations 

In the case of a single fluid with no entropy perturbations the pressure per- 
turbation 6p is simply a function of the density perturbation, 

5p = cl5p, (2.164) 

where the adiabatic sound speed is defined as 

c^ = ^. (2.165) 

In general the pressure perturbation (in any gauge) can be split into adia- 
batic and entropic (non-adiabatic) parts, by writing 

6p = cl5p + pT, (2.166) 

where F is the dimensionless entropy perturbation. The non-adiabatic part 
can be written as 

Sp^^,^pT^p'(^-^-^-^] . (2.167) 

The non-adiabatic pressure perturbation 5pnad, defined in this way is gauge- 
invariant, and represents the displacement between hypersurfaces of uniform 
pressure and uniform density. 

We can extend the notion of adiabatic and non-adiabatic from density and 
pressure perturbations to perturbations in general: An adiabatic perturbation 
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is one for which all perturbations 5x share a common value for 5x/x' , where 
x' is the time dependence of the background value of x. If the Universe is 
dominated by a single fluid with a definite equation of state, or by a single 
scalar field whose perturbations start in the vacuum state, then only adiabatic 
perturbations can be supported. If there is more than one fluid, then the 
adiabatic condition is a special case, but for instance is preserved if a single 
infiaton field subsequently decays into several components. However, pertur- 
bations in a second field, for instance the one into which the infiaton decays 
during preheating, typically violate the adiabatic condition. 

The entropy perturbation between any two quantities (which are spatially 
homogeneous in the background) has a naturally gauge-invariant definition, 
which follows from the obvious extension of Eq. (|2.167D , 



T,,^^A- (2-168) 

x' y' 

In a multi-component system the entropy perturbation F can be split into 
two parts 

r = ri,t + r,ei, (2.169) 

where Tint is the intrinsic and Frei the relative entropy perturbation. The 
pressure perturbation in the a fluid is given as 

5Pa = cl6pa + Pa^a , (2.170) 

where Fq, is intrinsic entropy perturbation of this particular fluid and the sound 
speed of the a fluid is 

cl = ^. (2.171) 
The sum of the intrinsic entropy perturbations in each fluid is represented by 

(2.172) 

a 

There may be some cases in which one can discuss the entropy perturbation 
within a single fluid, however we are more interested in the case where the 
entropy perturbation arises due to relative evolution between two or more 
fluids with different sound speeds. This is given by 

prrei = 5^(c2 -c^)5p«. (2.173) 

a 

The overall adiabatic sound speed is determined by the component parts[| 

= E# (2-174) 



^Note that there is an obvious typo on the right-hand-side of the first hne of Eq. (5.33) 
in Kodama and Sasaki [b5| . 
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Substituting this into equation (p.l73|) we obtain^ 
pTrei = \ E (^^i±M^^l±M(e2 _ cl)S^^ + + , (2.176) 



where, following Kodama and Sasaki we introduce new dimensionless 
variables for the perturbed total and relative energy densitie^ 

A = (2.177) 

P + P 

A„ ^ (2.178) 

Pa +Pa 

S^p = A^-Ap. (2.179) 

While this is a useful definition of the entropy perturbation Sa/s in the 
absence of energy transfer between the fluids (g^ = 0) it is not very well suited 
to (and not gauge invariant in) the more general case, as can be seen by the 
presence of adiabatic perturbation A in the expression for Frei- Instead it is 
useful to introduce alternative definitions 

A„ = 7^ TT^^^' (2-180) 

(1 - qa){pa+Pa) 

S^p = A^-Ap, (2.181) 

Saf3 is then a gauge invariant definition of the entropy perturbation, and it 
allows us to re-write Eq. ( p.l76| ) as 



ga)(Pa +Pa)(l - g/3)(P/3 + P/b) 

p + p 

which vanishes in the absence of any non-adiabatic perturbation 



2 P + P 

a,l3 F I 



2.11 Refined equations of motion 



Here we present the equations of motion in terms of Kodama and Sasaki's 
redefined density perturbations defined in our Eqs. (|2.177| - p.l79| ). After some 
delicate calculations we obtain from the equations presented in Section |2.10.5 
a system of coupled first order ordinary differential equations 



A' - 3/ic?A 



a — 



2/in ■ 
p + p. 



(2.183) 



cr 



ha 



A'KGa^{p + p) 
P -3k 



A + 



P 



P + P 



^There is a typo in Eq. (5.36) of Ref. The second term on the right-hand-side in 

their final equation for pFj-ei should be multiplied by A (in their notation). This agrees with 
the corrected version of the equations given in Ref. . 

^Note that Kodama and Sasaki use A to describe 5p/ p. 
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af3 



Va(3 + hv. 



where 



a/3 



+ 



n, 



3 p + p 

5p-2{k^ -3k)II/3' 



(2.184) 



-3h{qa - qp) 



P + P 



+2,h{E^p-T^p), (2.185) 
+hfa/3 - {clA^ - 4Afs) - - , (2.186) 



n 



a/3 



0/3 
0/3 



Va-V/B, 

Pa +Pa P/3 + P/3 

2{3k - F) / n„ 



Pa +Pa P/3+ Pl3 



Qa^a — Qp^P 
fa - fp- 



(2.187) 

(2.188) 

(2.189) 

(2.190) 
(2.191) 



The residual dependence on the individual and Va on the right-hand-sides 
of the equations of motion for S^p and v^p can be eliminated by substituting 
in 



p + p 



s, 



a7 ; 



E 



P7 +^7 
— 

P+P 



(2.192) 
(2.193) 



Finally then we obtain ^ 



S'ap +-h[qa{l + cl) + qp{l + cl)]Sap 

+ h[q^{l + cl)-qp{l + 4)]Y, 



P7 +P7 
P + P 



k'^Vap - Sh [c1{qa - qp) + ga(l + 4) - qp{l + 4)] A 
pT - 2{k^ - 3K)n/3 



V'ap +h 



-3h{qa - qp) 
■ 3 



P + P 



+ 3h{E^p-T^p), (2.194) 



t^a/3 



^These should be compared with Kodama and Sasaki's equations (5.53) and (5.57). There 
are some, possibly typographical, errors in their equations. Again our corrected equations 



agree with Ref. 
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-Ih [cl - 4 - + cl) + q^il + 4)] J2 ^II^^ Kt + 



= - c: 



4)A--(4 + 4)5«^ 



P7 +P7 



(S'a^ + Sp^) + hfalS - r„/3 - Ual3 .(2.195) 



These equations provide a closed set of coupled first-order equations in the 
absence of intrinsic entropy perturbation (F^ = 0) in the individual fluids, 
and if one specifies the energy and momentum transfer between fluids, and 
their anisotropic stress. 

However they remain rather unsatisfactory due to the presence of the adia- 
batic perturbation A as a source term for gf^ 7^ on the right-hand-side of the 
equation of motion for the entropy perturbation Sa/3- It should be possible to 
conceive of an adiabatic perturbation on large scales, Sf] = Spa/ p'„ for all fluids 
a (i.e. a perturbation along the homogeneous background trajectory), which 
remains adiabatic (i.e. remains along the trajectory) and does not generate 
an entropy perturbation even in the presence of energy transfer between the 
fluids. This should indeed be possible if we work in terms of our alternative 
(gauge-invariant for 7^ 0) entropy perturbation Sa/s- 

We can then write the equations of motion for Sa/3 and Va/s as 



al3 



1 
'2 
1 

+ 2 



+ 



l-qa i-qp 



- 3hqa{l + 4) - 3hq/3{l + c 



S, 



aj3 



l-qa 1 - qf3 

(1 - qy){pT+p^) 



-3%,(l + 4) + 3%(l + 4) 



7 

a/3 

3 



h 



P + P 



Sa^ + 3h{E. 



aj3 



1 - (^a + 4) + o?«(i + 4) + + 4) 



'Jap 



(2.196) 



^ ^, p ' p 



= -[4(1 - qa) - 4(1 - qp)]A - -[4(1 - qa) + 4(1 - qp)]Saf3 



4(1 



+hfa6 



- Qa) 
^al3 ~ i 



(1-?7)(P7+P 



P + P 



(>S'a7 + Sp^) 



(2.197) 



where 



Sa, = A-3V' 
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1 + 



UrrGa^p + p) 



Vaf3 
Eafi 



k"^ — 3/€ 



A + 3ha , 



1 - go 1 - g/3 

{1 - qa){pa + Pa) {I - qp){pp + Pp) 



(2.198) 
(2.199) 
(2.200) 
(2.201) 



and we have introduced the non-adiabatic part of the perturbed energy transfer 
for each fluid 

where Qa = 'ihqa{pa +Pa)/«- This vanishes if the perturbed energy transfer 
vanishes on surfaces of constant total density, and thus is zero if the energy 
transfer is a function solely of the total density, or if the perturbations are 
purely adiabatic. 

We have had to introduce the total density perturbation on uniform curva- 
ture hypersurfaces S'ak and the momentum perturbation on uniform curvature 
hypersurfaces Va — iSj/h). 

But there remains a source term on the right hand side of the S equation 
from adiabatic perturbations, S'ak, for spatially curved background FRW mod- 
els (k 7^ 0), which threatens to source entropy perturbations on large scales 
by adiabatic ones, in contradiction with our intuition. 



2.12 Conserved quantities on large scales 

During the different epochs of the universe different forms of matter dominate 
and may later become insignificant or decay altogether. Instead of following 
the evolution of the various matter fields it is more convenient to follow the 
evolution of a metric variable, the curvature perturbation, which is, as we shall 
see below, conserved on large scales for adiabatic perturbations on uniform 
density hypersurfaces. 

Bardeen, Steinhardt and Turner constructed a conserved quantity Cbst 
from quantities in the uniform expansion gauge, where the expansion scalar 
59 = 0, as defined in Eq. ( |2.31D , which can be identified with the curvature 
perturbation on uniform density hypersurfaces ( in 0, 0. A neat way of 
showing the constancy of ( is by using the energy conservation equation in the 
uniform density gauge, Eq. (|2.101|) . This approach, which does not depend on 
the Einstein field equations, was introduced first in [|105|| . 



As shown in Section p.l0.6| we can split the pressure perturbation into an 



adiabatic and a non-adiabatic part, 

5p = cl6p + 5p^^ , (2.203) 
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where 6pnad is the non-adiabatic pressure perturbation defined in Eq. ( p.l671 ). 
Note that the first term in the above equation is zero in the uniform density 
gauge and we can identify 5pnad with the pressure perturbation in the uniform 
density gauge, 6p5p. Rewriting Eq. ( p.l01| ) and using the fact that in the 
uniform density gauge ipSp = — C we find 

C = —5pn.a -V^v + E') . (2.204) 

P + P 



On sufficiently large scales gradient terms can be neglected and W^, ES 



C = —5pn^ . (2.205) 

p + p 

It follows that C is conserved on large scales for adiabatic perturbations, for 
which (5pnad = 0. We emphasize, that this result has been derived indepen- 
dently from the gravitational field equations, using only the conservation of 
energy |p.05|| ! Hence this result holds in any theory of relativistic gravity 



In Ref. the gauge-invariant variable ^mfb is defined as 

2 $' + $/i , , 

Cmpb = $ + TTTT— 2.206 
3(1 + w)h 

where w = po/po- On large scales (where we neglect spatial derivatives) and in 
fiat-space (k = 0) with vanishing anisotropic stresses (vr^^ = 0, which requires 
that $ = \1/ m^) all three quantities ip^n, V'p and Cmfb coincide. The curvature 
perturbation, in one or the other of these forms, is often used to predict the 
amplitude of perturbations re-entering the horizon scale during the radiation 
or matter dominated eras in terms of perturbations that left the horizon during 
an infiationary epoch, because they remain constant on super-horizon scales 
(whose comoving wavenumber k <^h) for adiabatic perturbations |^ . 



It is illustrative to investigate the coincidence of the three quantities ^/'m, 
Tpp and Cmfb niore closely. 

2.12.1 Single-field inflation 

The specific relation between the infiaton field and curvature perturbations 
depends on the choice of gauge. In practice the infiaton field perturbation 
spectrum can be calculated on uniform-curvature {ip = 0) slices, where the 
field perturbations have the gauge-invariant definition Q 

Sif^ = 5ifi + ^ij. (2.207) 

In the slow-roll limit the amplitude of field fiuctuations at horizon crossing 
(A = H^^) is given by H/2'7i (see Section |3.5.2| ). Note that this is the am- 



plitude of the asymptotic solution on large scales. This result is independent 

^Note that related results have been obtained in particular non-Einstein gravity theories 
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of the geometry and holds for a massless scalar field in de Sitter spacetime 
independently of the gravitational field equations. 

The field fluctuation is then related to the curvature perturbation on co- 
moving hypersurfaces (on which the scalar field is uniform, dipra = 0) using 
Eq. (^, by 

7^ = V^n, = -5^^ . (2.208) 

<^ 

We will now demonstrate that for adiabatic perturbations we can identify 
the curvature perturbation on comoving hypersurfaces, TZ, with the curvature 
perturbation on uniform-density hypersurfaces, —C,. In an arbitrary gauge 
the density and pressure perturbations of a scalar field are given by (see 
Eqs. i W^ ) and (VtTQD ) 

5p = a-\^'S^' ~(j)(p'^) + V^S^, (2.209) 
Sp = a-'^{^'6^' -(j)ip'^) -V^6ip, (2.210) 

where = dV/dip. Thus we find 6p—6p = 2V^6ip. The entropy perturbation, 
Eq. ( p.l67| ), for a scalar field is given by 



and from the energy and momentum constraints. Equations ( p.89|) and (|2.90| ), 
for a scalar field we get 

V'^{'4j-ha)= AttG {<p'6y^' - <p''^(j) - {<p" - h<p')6<p) . (2.212) 

Hence the entropy perturbation for a scalar field vanishes if gradient terms can 
be neglected. Therefore on large scales the scalar field perturbations become 
adiabatic and then on uniform-density hypersurfaces both the density and 
pressure perturbation must vanish and thus so does the field perturbation 
Sifsp = for 7^ 0. Hence the uniform- density and comoving hypersurfaces 
coincide, and TZ and —( are identical, for adiabatic perturbations and, from 
Eq. ( p.205|) , are constant on large scales. 



2.12.2 Multi-component infiaton field 

During a period of inflation it is important to distinguish between "light" fields, 
whose effective mass is less than the Hubble parameter, and "heavy" fields 
whose mass is greater than the Hubble parameter. Long-wavelength (super- 
Hubble scale) perturbations of heavy fields are under-damped and oscillate 
with rapidly decaying amplitude ((0^) oc a~^) about their vacuum expectation 
value as the universe expands. Light fields, on the other hand, are over-damped 
and may decay only slowly towards the minimum of their effective potential. It 
is the slow-rolling of these light fields that controls the cosmological dynamics 
during infiation. 
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The inflaton, defined as the direction of the classical evolution, is one of the 
light fields, while the other light fields (if any) will be taken to be orthogonal 
to it in field space. In a mult i- component inflation model there is a family of 
inflaton trajectories, and the effect of the orthogonal perturbations is to shift 
the inflaton from one trajectory to another. 

If all the fields orthogonal to the inflaton are heavy then there is in effect a 
unique inflaton trajectory in field space. In this case even a curved path in field 
space, after canonically normalizing the inflaton trajectory, is indistinguishable 
from the case of a straight trajectory, and leads to no variation in (. 

When there are multiple light fields evolving during inflation, uncorrelated 
perturbations in more than one field will lead to different regions that are 
not simply time translations of each other. In order to specify the evolution 
of each locally homogeneous universe one needs as initial data the value of 
every cosmologically significant field. In general, therefore, there will be non- 
adiabatic perturbations, T^y 7^ 0. 

If the local integrated expansion, A^, is sensitive to the value of more than 
one of the light fields then ( is able to evolve on super-horizon scales, as 



has been shown by several authors |9^, Note also that the comoving 



and uniform- density hypersurfaces need no longer coincide in the presence of 
non-adiabatic pressure perturbations. In practice it is necessary to follow the 
evolution of the perturbations on super-horizon scales in order to calculate 
the curvature perturbation at later times. In most models studied so far, the 
trajectories converge to a unique one before the end of inflation, but that need 
not be the case in general. 

The separate universe approach described in Section |2.13| gives a rather 



straightforward procedure for calculating the evolution of the curvature per- 
turbation, ip, on large scales based on the change in the integrated expansion, 
N, in different locally homogeneous regions of the universe. This approach 
was developed in Refs. fO^, for general relativistic models where scalar 
fields dominate the energy density and pressure, though it has not been ap- 
plied to many specific models. In the case of a single-component inflaton, 
this means that on each comoving scale. A, the curvature perturbation, (, on 
uniform- density (or comoving orthogonal) hypersurfaces must stop changing 
when gradient terms can be neglected. More generally, with a mult i- component 
inflaton, the perturbations generated in the fields during inflation will still de- 
termine the curvature perturbation, (, on large scales, but one needs to follow 
the time evolution during the entire period a scale remains outside the horizon 
in order to evaluate ( at later times. This will certainly require knowledge of 
the gravitational field equations and may also involve the use of approxima- 
tions such as the slow-roll approximation to obtain analytic results. 

2.12.3 Preheating 

During inflation, every fleld is supposed to be in the vacuum state well before 
horizon exit, corresponding to the absence of particles. The vacuum fluctu- 
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ation cannot play a role in cosmology unless it is converted into a classical 
perturbation, defined as a quantity which can have a well-defined value on a 



sufficiently long time-scale 71 1. For every light field this conversion occurs 
at horizon exit (A ~ H~^). In contrast, heavy fields become classical, if at all, 
only when their quantum fiuctuation is amplified by some other mechanism. 

There has recently been great interest in models where vacuum fiuctuations 
become classical (i.e., particle production occurs) due to the rapid change 
in the effective mass (and hence the vacuum state) of one or more ffelds. 
This usually (though not always ||T^) occurs at the end of inffation when 
the inffaton oscillates about its vacuum expectation value which can lead to 
parametric amplification of the perturbations — a process which has become 
known as preheating |5^. The rate of amplification tends to be greatest for 
long-wavelength modes and this has lead to the claim that rapid amplification 
of non-adiabatic perturbations could change the curvature perturbation, ^, 
even on very large scales Preheating will be discussed in more detail in 

Section ^. 

Within the separate universes picture, discussed in Section p.l3| , this is cer- 
tainly possible if preheating leads to different integrated expansion in different 
regions of the universe. In particular by Eq. ( p.205| ) the curvature perturbation 
C can evolve if a significant non-adiabatic pressure perturbation is produced 
on large scales. However it is also apparent in the separate universes picture 
that no non-adiabatic perturbation can subsequently be introduced on large 
scales if the original perturbations were purely adiabatic. This is of course also 
apparent in the field equations where preheating can only amplify pre-existing 
field fiuctuations. 

Efficient preheating requires strong coupling between the infiaton and pre- 
heating fields which typically leads to the preheating field being heavy during 
infiation (when the infiaton field is large). The strong suppression of super- 
horizon scale fiuctuations in heavy fields during infiation means that in this 
case no significant change in C, is produced on super-horizon scales before 
back-reaction due to particle production on much smaller scales damps the 
oscillation of the infiaton and brings preheating to an end [^, |62[ . 

Because the first-order effect is so strongly suppressed in such models, 
the dominant effect actually comes from second-order perturbations in the 
fields 49, 62]. The expansion on large scales is no longer independent of 
shorter wavelength field perturbations when we consider higher-order terms in 
the equations of motion. Nonetheless in many cases it is still possible to use 
linear perturbation theory for the metric perturbations while including terms 
quadratic in the matter field perturbations. In Ref. this was done to show 
that even allowing for second-order field perturbations, there is no significant 
non-adiabatic pressure perturbation, and hence no change in on large scales 
in the original model of preheating in chaotic infiation. 

More recently a modified version of preheating has been proposed [0] (re- 
quiring a different model of inflation) where the preheating field is light during 
inflation, and the coupling to the infiaton only becomes strong at the end of 
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inflation. In such a multi-component inflation model non-adiabatic perturba- 
tions are no longer suppressed on super-horizon scales and it is possible for 
the curvature perturbation ( to evolve both during inflation and preheating, 
as described in Section p.l2.2| . 

2.12 A Non-interacting multi-fluid systems 

In a multi-fluid system we can define uniform-density hypersurfaces for each 
fluid and a corresponding curvature perturbation on these hypersurfaces, 

Cm^-^-^/^- (2.213) 

Equation ( p.204| ) then shows that remains constant for adiabatic pertur- 
bations in any fluid whose energy-momentum is locally conserved: 

^^^4 = • (2.214) 

Using the relation between the curvature perturbation on constant density hy- 
persurfaces and the density perturbation on constant curvature hypersurfaces 

h ~ 

^<Sp = —,^P^ , (2.215) 
p' 

the total curvature perturbation, i.e. curvature perturbation on uniform total 
density hypersurfaces, can then be defined as 

C=-^^^^, (2.216) 

2^ a Pa 

where the sum is over the fluid components. This quantity is not in general 
conserved in a multi-fluid system in presence of a non-zero relative entropy 
perturbation Frei 7^ in Eq. ( ^.1821 ). 



Thus, for example, in a universe containing non-interacting cold dark mat- 
ter plus radiation, which both have well-defined equations of state (pm = and 

= P7/3), the curvatures of uniform- matter-density hypersurfaces. Cm, and of 
uniform-radiation-density hypersurfaces, remain constant on super- horizon 
scales. The curvature perturbation on the uniform-total-density hypersurfaces 
is given by 

= (4/3)p,C. + PmCm _ 

(4/3)p, + p„, ^ ^ 

At early times in the radiation dominated era (p^ ^ pm) we have Cini — C7; 
while at late times (pm S> p-y) we have Cfin — Cm. C remains constant through- 
out only for adiabatic perturbations where the uniform-matter-density and 
uniform-radiation-density hypersurfaces coincide, ensuring C7 = Cm. The 
isocurvature (or entropy) perturbation is conventionally denoted by the per- 
turbation in the ratio of the photon and matter number densities 

Sn 

S-ym = = 3 (C7 — Cm) . (2.218) 
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Hence the entropy perturbation for any two non-interacting fluids always re- 
mains constant on large scales independent of the gravitational field equations. 
Hence we recover the standard result for the final curvature perturbation in 
terms of the initial curvature and entropy perturbation 



,fin Sini 




(2.219) 



2.13 The separate universe picture 

Thus far we have shown how one can use the perturbed field equations, to 
follow the evolution of linear perturbations in the metric and matter fields 
in whatever gauge one chooses. This allows one to calculate the correspond- 
ing perturbations in the density and pressure and the non-adiabatic pressure 
perturbation if there is one, and see whether it causes a significant change in 

C- 

However, there is a particularly simple alternative approach to studying 
the evolution of perturbations on large scales, which has been employed in 
some multi-component inflation models [^, ^ |9^, 98, dl, 75|. This considers 



each super-horizon sized region of the Universe to be evolving like a separate 
Robert son- Walker universe where density and pressure may take different val- 
ues, but are locally homogeneous. After patching together the different re- 
gions, this can be used to follow the evolution of the curvature perturbation 
with time. Figure 1 shows the general idea of the separate universe picture, 
though really every point is viewed as having its own Robert son- Walker region 
surrounding it. 

Consider two such locally homogeneous regions (a) and (6) at fixed spatial 
coordinates, separated by a coordinate distance A, on an initial hypersurface 
(e.g., uniform- density hypersurface) specified by a fixed coordinate time, t = 
ti, in the appropriate gauge (e.g., uniform-density gauge). The initial large- 
scale curvature perturbation on the scale A can then be defined (independently 
of the background) as 

6^P^ = ^Pai - ^61 . (2.220) 

On a subsequent hypersurface defined by t = ^2 the curvature perturbation at 
(a) or (6) can be evaluated using Eq. ( p.45|) [but neglecting V^a] to give 



^a2 = i^al-5Na, (2.221) 

where the integrated expansion between the two hypersurfaces along the world- 
line followed by region (a) is given by A^^a = + ^A^, and A^ = Ina is the 
expansion in the unperturbed background and 

/■'I ~ 

6K = j -dOadt. (2.222) 
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This result was derived first by solving a differential equation |p6[ , and then [TTa by inte- 



grating Eq. (2.204) using Eq. (2.217). We have here demonstrated that even the integration 



is unnecessary. 
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The curvature perturbation when t = t2 on the comoving scale A is thus given 
by 

S^p2 = ^a2 - = 5^1 - {Na - iVfe) . (2.223) 

In order to calculate the change in the curvature perturbation in any gauge on 
very large scales it is thus sufficient to evaluate the difference in the integrated 
expansion between the initial and final hypersurface along different world-lines. 
In particular, using Eq. (|2.223| ), one can evolve the curvature perturba- 



tion, on super-horizon scales, knowing only the evolution of the family of 
Robertson-Walker universes, which according to the separate Universe as- 
sumption describe the evolution of the Universe on super-horizon scales: 

-AC = -AAr. (2.224) 

where —C, = ipa — 4'b on uniform- density hypersurfaces and AA^ = Na — Nb in 
Eq. ( |2.223| ). As we shall discuss in Section ^, this evolution is in turn specified 
by the values of the relevant fields during inflation, and as a result one can 
calculate ( at horizon re-entry from the vacuum fluctuations of these fields. 

While it is a non-trivial assumption to suppose that every comoving re- 
gion well outside the horizon evolves like an unperturbed universe, there has 
to be some scale As for which that assumption is true to useful accuracy. 
If there were not, the concept of an unperturbed (Robertson- Walker) back- 
ground would make no sense. We use the phrase 'background' to describe 
the evolution on a much larger scale Aq, which should be much bigger even 
than our present horizon size, with respect to which the perturbations in this 
section are defined. It is important to distinguish this from regions of size As 
large enough to be treated as locally homogeneous, but which when pieced 
together over a larger scale. A, represent the long- wavelength perturbations 
under consideration. Thus we require a hierarchy of scales: 

Ao > A > As > cH-^ . (2.225) 

Ideally Aq would be taken to be infinite. However it may be that the Universe 
becomes highly inhomogeneous on some very much larger scale, Ae ^ Aq, 
where effects such as stochastic or eternal inflation determine the dynamical 
evolution. Nevertheless, this will not prevent us from defining an effectively 
homogeneous background in our observable Universe, which is governed by the 
local Einstein equations and hence impervious to anything happening on vast 
scales. Specifically we will assume that it is possible to foliate spacetime on 
this large scale Aq with spatial hypersurfaces. 

When we use homogeneous equations to describe separate regions on length 
scales greater than As, we are implicitly assuming that the evolution on these 
scales is independent of shorter wavelength perturbations. This is true within 
linear perturbation theory in which the evolution of each Fourier mode can 
be considered independently, but any non-linear interaction introduces mode- 
mode coupling which undermines the separate universes picture. The separate 
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universe model may still be used for the evolution of linear metric perturbation 
if the perturbations in the total density and pressure remain small, but a 
suitable model (possibly a thermodynamic description) of the effect of the 
non-linear evolution of matter fields on smaller scales may be necessary in 
some cases. 

Adiabatic perturbations in the density and pressure correspond to shifts 
forwards or backwards in time along the background solution, 5p/5p = p' /p' = 
cj., and hence F = in Eq. ( p. 167 ). For example, in a universe containing only 
baryonic matter plus radiation, the density of baryons or photons may vary 
locally, but the perturbations are adiabatic if the ratio of photons to baryons 
remains unperturbed (cf. Eq.( p.21^ )). Different regions are compelled to un- 
dergo the same evolution along a unique trajectory in field space, separated 
only by a shift in the expansion. The pressure p thus remains a unique function 
of the density p and the energy conservation equation, dp/dN = —3{p + p), 
determines p as a function of the integrated expansion, N. Under these con- 
ditions, uniform-density hypersurfaces are separated by a uniform expansion 
and hence the curvature perturbation, (, remains constant. 

For F 7^ it is no longer possible to define a simple shift to describe both 
the density and pressure perturbation. The existence of a non-zero pressure 
perturbation on uniform-density hypersurfaces changes the equation of state 
in different regions of the Universe and hence leads to perturbations in the 
expansion along different worldlines between uniform-density hypersurfaces. 
This is consistent with Eq. ( |2.205D which quantifies how the non-adiabatic 



pressure perturbation determines the variation of ( on large scales [75[] . 

We defined a generalized adiabatic condition which requires Txy = for 
any physical scalars x and y in Eq. (|2.168|) in Section p.l0.6| . In the separate 
universes picture this condition ensures that if all field perturbations are adia- 
batic at any one time (i.e. on any spatial hypersurface), then they must remain 
so on large scales at any subsequent time. Purely adiabatic perturbations can 
never give rise to entropy perturbations on large scales as all fields share the 
same time shift, 6r] = 6x/x', along a single phase-space trajectory. 
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3 Quantum fluctuations from inflation 



In this section we will discuss the generation of quantum fluctuations from 
inflation in the context of a particularly simple mult i- component model: as- 
sisted inflation. The assisted inflation model is one of the very few models, in 
which it is possible to calculate the power spectra of the homogeneous and the 
inhomogeneous field fiuctuations analytically, without resorting to the slow 
roll approximation |Bo|, |75|, |61|. 



3.1 The assisted inflation model 

A single scalar field with an exponential potential is known to drive power-law 
infiation, where the cosmological scale factor grows as a oc with p > 1, for 
sufficiently fiat potentials |]69| , |14| , p!9| . Liddle, Mazumdar and Schunck |^ 



recently proposed a novel model of infiation driven by several scalar fields with 
exponential potentials. Although each separate potential. 



Vi = Vq exp 



Pi mpi 



(3.1) 



may be too steep to drive infiation by itself {pi < 1), the combined effect of 
several such fields, with total potential energy 



i=l 



leads to a power-law expansion a oc with |^| 



(3.2) 



P = ^Pi, 

1=1 



(3.3) 



provided p > 1/3. Supergravity theories typically predict steep exponential 
potentials, but if many fields can cooperate to drive infiation, this may open 
up the possibility of obtaining infiationary solutions in such models. 

Scalar fields with exponential potentials are known to possess self-similar 
solutions in Friedmann- Robert son- Walker models either in vacuum H\ or 
in the presence of a barotropic fiuid ||106| , |103| , p!9| , In the presence of 



other matter, the scalar field is subject to additional friction, due to the larger 
expansion rate relative to the vacuum case. This means that a scalar field, 
even if it has a steep (non-infiationary) potential may still have an observable 



100 



dynamical effect in a radiation or matter dominated era [ 107 , 25, 2C, 

The paper of Liddle, Mazumdar and Schunck [^] was the first to consider 
the effect of additional scalar fields with independent exponential potentials. 
They considered n scalar fields in a spatially fiat Friedmann-Robertson- Walker 
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universe with scale factor a{t). The Lagrange density for the fields, as an 
extension of Eq. (|2.113|) , is 



n 



(3.4) 



i=l 



with each exponential potential Vi of the form given in Eq. ( |3.1[ ). The cos- 
mological expansion rate is then given by the Friedmann equation (|2.85|) in 
coordinate time 

^ 8- - (3.5) 



3mpj 



|:(..i.?) 



and the individual fields obey the Klein-Gordon equation (|2.120|) 



(3.6) 



One can then obtain a scaling solution of the form [32 



(3.7) 



Differentiating this expression with respect to time, and using the form of the 
potential given in Eq. (p.l|) then implies that 



1 



and hence 



The scaling solution is thus given by |63| 



mpi 



Vie 



TC 



In^. 

Pi 



(3.8) 



(3.9) 



(3.10) 



A numerical solution with four fields is shown in Fig. ^ as an example. In 
Ref. [33 1 the authors demonstrated the existence of a scaling solution for n 
scalar fields written in terms of a single re-scaled field (p = \/p/pi^i- The 
choice of ipi rather than any of the other fields is arbitrary as along the scaling 
solution all the (pi fields are proportional to one another. 

In this part of the thesis we will prove that this scaling solution is the late- 
time attractor by choosing a redefinition of fields (a rotation in field space) 
which allows us to write down the effective potential for field variations or- 
thogonal to the scaling solution and show that this potential has a global 
minimum along the attractor solution. In general the full expression for an 
arbitrary number of fields is rather messy so we first give, in Section |3.2|, the 



simplest case where there are just two fields, and then extend this to n fields 



40 



0.1 



10 
time, t 



100 



1000 



Figure 2: Evolution of four fields ipi, ip2, v^s and ip4 (from bottom to top) 
during assisted inflation with pi = 0.3, p2 = 1, Pa = 2 and p4 = 7. 



in Section p.3| . The resulting inflationary potential is similar to that used in 
models of hybrid inflation and we show in Section p.4| that assisted inflation 
can be interpreted as a form of "hybrid power-law inflation" . As in the case of 
power-law or hybrid inflation, one can obtain analytic expressions for inhomo- 
geneous linear perturbations close to the attractor trajectory without resorting 
to slow-roll type approximations. Thus we are able to give exact results for 
the large-scale perturbation spectra due to vacuum fluctuations in the fields 
in Section We discuss our results in Sections ^]6| and |^. 



3.2 Two field model 

We will restrict our analysis initially to just two scalar fields, (pi and ip2, with 
the Lagrange density 



exp 



We define the fields 

^2 = 



+ 



VP1+P2 

mpi 



0"2 



i/l67r(pi +P2) 



Pi rripi 



exp 



P2 rnpi 



(3.11) 



Pi In 



Pi 



Pi +P2 



+ P2 In 



P2 



Pi +P2 



VPi + P2 



+ 



mpi 



In: 



vT67r y P1+P2 P2 



(3.12) 
(3.13) 
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Figure 3: The field rotation in the two field case. The ^2 direction is along 
the attractor, orthogonal to it. 



to describe the evolution along and orthogonal to the scaling solution, re- 
spectively, by applying a Gram-Schmidt orthogonalisation procedure. The 
re-defined fields (p2 and are orthonormal linear combinations of the original 
fields and ^f^- They represent a rotation, and arbitrary shift of the origin, 
in field-space. Fig. |] illustrates this. Thus <^2 and have canonical kinetic 
terms, and the Lagrange density given in Eq. ( |3.11| ) can be written as 



i(Va2')-^(t^2)exp 



pi + p2 mpi 



(3.14) 



where 



V{a2] 



Vn 



+ 



Pi 



Pi + P2 
P2 



exp 



Pi +P2 



exp 




(3.15) 



It is easy to confirm that V{a2) has a global minimum value Vq at ct2 = 0, 
which implies that ct2 = is the late time attractor, which coincides with the 
scaling solution given in Eq. ( glOD for two fields. 

Close to the scaling solution we can expand about the minimum, to second- 
order in o"2, and we obtain 



V{^2,(T2) ~ Vo 



1 + 



Stt 



- 2 
0"o 



(Pi +P2) ml^. 



exp 



IGtt 



V52 



Pi + p2 mpi 



(3.16) 
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Note that the potential for the field a2 has the same form as in models of 
hybrid infiation ||6^, ^ where the infiaton field rolls towards the minimum 
of a potential with non-vanishing potential energy density Vq. Here there is in 
addition a "dilaton" field, (f2, which leads to a time-dependent potential energy 
density as o"2 — >■ 0. Assisted infiation is related to hybrid infiation p6|, |67|, pO 



in the same way that extended infiation [£8| was related to Guth's old infiation 
model 



m 



As in hybrid or extended infiation, we require a phase transition 
to bring infiation to an end. Otherwise the potential given by Eq. ( p. 16 ) leads 
to infiation into the indefinite future. 



3.3 Many field model 

We will now prove that the attractor solution presented in Ref. |^ is the global 
attractor for an arbitrary number of fields with exponential potentials of the 
form given in Eq. (|3.1| ), using proof by induction. To do this, we recursively 
construct the orthonormal fields and their potential. 

Let us assume that we already have n fields (fi with exponential potentials 
Vi of the form given in Eq. ( ^.l] ) and that it is possible to pick n orthonormal 
fields a2, ■ ■ ■ ,(Tn and such that the sum of the individual potentials Vi can 
be written as 

±V, = V^e.J-.l^^), (3.17) 

where we will further assume that Vn = Ki(a'i) has a global minimum \4(0) = 
Vq when cTj = for all i from 2 to n. 

It is possible to extend this form of the potential to n+1 fields if we consider 
an additional field fn+i with an exponential potential Ki+i of the form given 
in Eq. ( p.l|) . Analogously to the two field case, we define 

VKVn + y/Pn+lVn+1 



VPn+1 
?1 

'I67r(p„+i) V " Pn+l ' ' Pn+1 



+ , ^ =hPnln- \-pn+i\n- , (3.18) 



where 



^/p^^n - VK^n+l , mpi PnPn+1 , Pn 

an+1 = 1== + ^=W— In , 3.19 

VPn+1 VlQn V Pn+1 Pn+1 



Pn+1 =Pn+ Pn+1 ■ (3.20) 



Using these definitions we can show that the sum of the n + 1 individual 
potentials Vi can be written as 



n+1 



i=l 




j:V = Vn,,e.pi-J^^] , (3.21) 



43 



where K+i = 14+i(crj) is given by 

Pr, 



Vr 



n+1 



Vr 



n — 
Pn+1 



exp 



+ Vo exp 

Pn+1 




(3.22) 



Because we have assumed that Vn has a global minimum value V^(0) = Vq 
when cjj = for all i from 2 to n, one can verify that Vn+i also has a minimum 
value Ki+i(0) = Vo when o"j = 0, for all i from 2 to n + 1. 

However, we have already shown in Section |3]^ that for two fields and 
ip2, we can define two fields (^2 and o"2, given in Eqs. ( p.l2| ) and ( p.l3| ) whose 
combined potential given in Eq. ( 3.14 ) is of the form required in Eq. ( 3.17| ), 
with P2 = Pi + P2- Hence we can write the potential in the form given in 
Eq. ( p.l7|) for n fields, for all n > 2, with 



P=Pn 



Y^Pi 



(3.23) 



i=l 



Equations ( p.l2| ) and ( |3.18| ) then lead us to the non-recursive expression for 
the "weighted mean field" 



E 

i=l 



Pi 



, Pi 
~ 

/Ibirp p 



(3.24) 



which describes the evolution along the scaling solution. This is simply a 
rotation in field space plus an arbitrary shift, chosen to preserve the form of 
the potential given in Eq. (|3.17| ). The n — 1 fields cxj describe the evolution 
orthogonal to the attractor trajectory. 

The potential Vn has a global minimum at (Xj = 0, which demonstrates 
that this is the stable late-time attractor. From Eqs. ( ^.ISp and ( |3.22| ) we get 
a closed expression for Vn, 



Vn 



+ 




(3.25) 



Close to the attractor trajectory (to second order in cij) we can write a Taylor 
expansion for the potential 




exp 



J=2 



167r (fin 

p rripi 



(3.26) 
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Note that this expression is dependent only upon p and not on the individual 
Pi- 



3.4 Stringy hybrid inflation 



The form of the potentials in Eqs. (|3.16 ) and ( |3.26|) is reminiscent of the 
effective potential obtained in the Einstein conformal frame from Brans-Dicke 
type gravity theories ||102|| . The appearance of the weighted mean field, as 
a "dilaton" field in the potential suggests that the matter Lagrangian might 
have a simpler form in a conformally related frame. If we work in terms of a 
conformally re-scaled metric 



gf,^ = exp 



p mpi 



9^J.l' ■ 



(3.27) 



then the Lagrange density given in Eq. ( |3.14| ) becomes 



exp 




(3.28) 



In this conformal related frame the field (p is non-minimally coupled to 
the gravitational part of the Lagrangian. The original field equations were de- 
rived from the full action, including the Einstein-Hilbert Lagrangian of general 
relativity. 



S 



d Xy 



-9 



m 



PI 



16n 



R + C 



(3.29) 



where R is the Ricci scalar curvature of the metric g^j,^. In terms of the 
conformally related metric given in Eq. ( |3.27D this action becomes (up to 



boundary terms ||104]| ) 



S 



m 



PI 



IGtt 



^ n—l 

R-UJ{VfAil)--J2i^^^f-V 



i=l 



where we have introduced the dimensionless dilaton field 



IGvr if 



p mpi 

and the dimensionless Brans-Dicke parameter 

p — 3 



UJ 



(3.30) 



(3.31) 



(3.32) 



Thus the assisted inflation model is identical to — 1 scalar fields o"j with 
a hybrid inflation type potential V{ai) in a string-type gravity theory with 
dilaton, yj^ii oc ip. However, we note that in order to obtain power-law inflation 
with p ^ 1 the dimensionless constant u must be much larger than that found 
in the low-energy limit of string theory where u = —1. 
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3.5 Perturbations about the attractor 



The redefined orthonormal fields and the potential allow us to give the equa- 
tions of motion for the independent degrees of freedom. If we consider only 
linear perturbations about the attractor then the energy density is indepen- 
dent of all the fields except (p^ and we can solve the equation for the ^ field 
analytically. 

The field equation for the weighted mean field is 



167r V 



p mpi 

Along the line cij = for all i in field space, we have 



V = Vq exp 



IGtt (p 
p mpi 



(3.33) 



(3.34) 



tractor ||36|, for this potential, where 



and the well-known power-law solution ||69[ with a oc t*^ is the late-time at- 

(3.35) 



^0 



and (fo = mpiA/p/47r and = mpiA/p/87rVo(3j9 - 

3.5.1 Homogeneous linear perturbations 

The field equations for the fields are 



dV 

(Ji ■\- 3Hai + — 
oai 



0. 



(3.36) 



where the potential V is given by Eqs. ( p.l7| ) and (|3.25|) , and the attractor 
solution corresponds to cij = 0. Equation ( |3.26| ) shows that we can neglect 
the back-reaction of o"j upon the energy density, and hence the cosmological 
expansion, to first-order and the field equations have the solutions 



where 



3p-l 



-1 ± 



'3(p-3) 



3p 



(3.37) 



(3.38) 



for p > 3, confirming that = is indeed a local attractor. In the limit p ^ oo 
we obtain Si = —2. For 1 < p < 3 the perturbations are under-damped and 
execute decaying oscillations about o"j = 0. 

The form of the solutions given in Eq. ( |3.37| ) for a'j(t) close to the attractor 
is the same for all the orthonormal fields cxj, as demonstrated in Fig. ^ Their 
evolution is independent of the individual Pi and determined only by the sum, 
p, as expected from the form of the potential given in Eq. ( p.26| ) . 
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Figure 4: Evolution of the fields ai, (T2, and as, orthogonal to the scaling 
solution, in the assisted inflation model shown in Fig. 0. 



3.5.2 Inhomogeneous Linear Perturbations 

Conventional hybrid inflation and power-law inflation are two of the very few 
models pO] in which one can obtain exact analytic expressions for the spectra 



of vacuum fluctuations on all scales without resorting to a slow-roll type ap- 
proximation. In the case of hybrid inflation, this is only possible in the limit 
that the inflaton field a approaches the minimum of its potential and we can 
neglect its back- reaction on the metric |3^. As the present model is so closely 



related to both power-law and hybrid inflation models close to the attractor, 
it is maybe not surprising then that we can obtain exact expressions for the 
evolution of inhomogeneous linear perturbations close to the scaling solution. 

We will work in terms of the redefined fields and (Tj, and their pertur- 
bations on spatially flat hypersurfaces ||7^. In the limit that o"j ^ we can 



neglect the back-reaction of the Oi field upon the metric and the field t^. Per- 
turbations in the field then obey the usual equation for a single field driving 



inflation and perturbations in the field Gi evolve in a fixed background. 
Defining 

u = aSif) , (3.39) 
Vi = a5ai , (3.40) 

we obtain the decoupled equations of motion for perturbations with comoving 
wavenumber fc, 

ul+(e-^-\uk = 0, (3.41) 
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+ = 0, (3.42) 



where z = a^^p'/a' and a prime denotes differentiation with respect to 
conformal time rj = j dt/a. For power-law expansion we have 2; oc a oc 
(_^)-p/(p-i) and thus 

a" p(2p-l) . 

We also have aH oc —p/ {{p — 1)1]) which gives 

,d^V _ 2(3p-l) _2 



da"^ {p 



r\ (3.44) 



where we have used the fact that d'^V/daf = IGnV/mpi along the attractor. 
The equations of motion therefore become 

,2 



= 0, (3.45) 
+ { - ^^^^y.k = 0, (3.46) 



where 



A = I^^EMEM, (3.48) 

and the general solutions in terms of Hankel functions are 

Uk = U,{-kr]y/^Hi'\-kr]) + U^^-kr^f'^ H'^K-kT^) , (3.49) 

Vik = Vu{-kr^f'^H^l\-kr^) + V2^{-kr^f'^H'^^\-kr^) . (3.50) 

Taking only positive frequency modes in the initial vacuum state for \kri\ ^ 1 
and normalising requires Uk and Vik — > e"*'^''/ v^2fc, which gives the vacuum 
solutions 



-(-7rr/)i/2gf(.+i)i^(i)(_^^)^ ^35^) 



2 



i(_vrry)l/2ef(A+l)^^(l)(_^^)^ (3 52) 



v.k - 2 

In the opposite limit, i.e., —kr] 0, we use the limiting form of the Hankel 
functions, hIP{z) ~ — (i/7r)r(z/)z~'^, and therefore on large scales, and at late 
times, we obtain 

Uk '^e'^^'^i-kr^Y-'Viv), (3.53) 
^e^t"(-A;r7)^-^r(A). (3.54) 
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The power spectrum of a Gaussian random field X is conventionally given 
by 161|] 

(1^) \Xk\'. (3.55) 

The power spectra on large scales for the field perturbations Sep and Sai are 
thus 



C(z/) H_ 



V = 7^^i-kvr^^', (3.56) 



where we have used r] = —{u — l/2)/{aH) and we define 

C{a) = (3.58) 

Both the weighted mean field cp and the orthonormal fields d"j are "light" fields 
(m^ < 3if^/2) during assisted inflation (for p > 3) and thus we obtain a 
spectrum of fluctuations in all the fields on large scales. Note that in the de 
Sitter limit, p — > oo and thus — >• 3/2 and A 3/2, we have P^^^ — > H/2tt, 

and ///27r. 

Denoting the scale dependence of the perturbation spectra by 

An^ ^ , (3.59) 



we obtain 



Ans^ = 3-2u = , (3.60) 

p — 1 



A„,,, = 3-2A = 3h- V^(ZEMEM|. (3.61) 
3.6 Conclusions 

We have shown that the assisted inflation model, driven by many scalar fields 
with steep exponential potentials, can be better understood by performing a 
rotation in field space, which allows us to re- write the potential as a product of 
a single exponential potential for a weighted mean field, (f, and a potential Vn 
for the orthogonal degrees of freedom, cxj, which has a global minimum when 



cXi = 0. This proves that the scaling solution found in Ref. is indeed the 
late-time attractor. 

The particular form of the potential which we present for scalar fields 
minimally-coupled to the spacetime metric, can also be obtained via a confor- 



mal transformation of a hybrid inflation type inflationary potential [^, ^ ^| 
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with a non-minimally coupled, but otherwise massless, dilaton field, yjdii oc if. 
Thus we see that assisted infiation can be understood as a form of power- 
law hybrid infiation, where the false-vacuum energy density is diluted by the 
evolution of the dilaton field. 

We have also been able to give exact solutions for inhomogeneous linear 
perturbations about the attractor trajectory in terms of our rotated fields. 
Perturbations in the weighted mean field (f corresponds to the perturbations 
in the density on the uniform curvature hypersurfaces, or equivalently, pertur- 
bations in the curvature of uniform density hypersurfaces: 

( = -^S!f. (3.62) 

These perturbations are along the attractor trajectory, and hence describe 
adiabatic curvature perturbations. The perturbation spectrum for ( is 

Vc = ij] ^<5^ (3-63) 

which will be V^^"^ = H'^/{2TT(f) in the de Sitter limit (where p oo). The 
spectral index of the curvature perturbations on large scales is thus given from 
Eq. ( CT) as 



dlnVc 2 

and is always negatively tilted with respect to the Harrison-Zel'dovich spec- 
trum where = 1. Note that in the de Sitter limit (p — oo) we recover the 
result of Ref. . 

First-order perturbations in the fields orthogonal to the weighted mean field 
are isocurvature perturbations during infiation. Vacuum fiuctuations lead to 
a positively tilted spectrum. The presence of non-adiabatic perturbations can 
lead to a more complicated evolution of the large-scale curvature perturbation 



than may be assumed in single-field infiation models |92, 90, as discussed 



in Section |^ and in particular Section |2.12| . However, we have shown that 



these perturbations decay relative to the adiabatic perturbations and hence 
we recover the single field limit at late times. In particular we find that 
the curvature perturbation ( becomes constant on super-horizon scales during 
infiation. Note, however, that assisted infiation must be ended by a phase 
transition whose properties are not specified in the model. If this phase tran- 
sition is sensitive to the isocurvature (non-adiabatic) fiuctuations orthogonal 
to the attractor trajectory, then the curvature perturbation, (, during the sub- 
sequent radiation dominated era may not be simply related to the curvature 
perturbation during infiation. 

Since the publication of the first two papers on assisted infiation, [75| , 
the model has enjoyed considerable further interest by several authors. In 
Ref. Copeland, Mazumdar and Nunes extended the model presented 



above to include cross-couplings between the scalar fields. Kanti and Olive 
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proposed a realisation of assisted inflation based on the compactification of a 
five-dimensional Kaluza-Klein model, extending the model to include standard 
chaotic-type potentials Green and Lidsey showed that an arbitrary 



number of scalar fields with exponential potentials can be obtained from com- 
pactification of a higher dimensional theory |^ . They used a rotation in field 
space to show that the system reduces to a single scalar field with a single 
exponential potential. Coley and van den Hoogen investigated the dynam- 
ics of assisted inflation, including a barotropic fluid and curvature as well as 
scalar fields with exponential potentials, applying methods of dynamical sys- 



tems ||T8[. Lately Kaloper and Liddle used the assisted inflation mechanism 
in connection with chaotic inflation and showed that the dynamics does not 
become independent of the initial conditions at late times ||51|| . 
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4 Preheating 



The standard inflationary paradigm is an extremely successful model in ex- 
plaining observed structures in the Universe (see Refs. pO, lUTI, ^ for reviews). 



The inhomogeneities originate from the quantum fluctuations of the inflaton 
field, which on being stretched to large scales become classical perturbations 
as shown in the preceding section. The field inhomogeneities generate a per- 
turbation in the curvature of uniform density hypersurfaces, (, as shown in the 
preceding section and later on these inhomogeneities are inherited by matter 
and radiation when the inflaton field decays. In the simplest scenario, the cur- 
vature perturbation on scales much larger than the Hubble length is constant, 
and in particular is unchanged during the inflaton decay. This enables a pre- 
diction of the present-day perturbations which does not depend on the specific 
cosmological evolution between the late stages of inflation and the recent past 
(say, before nucleosynthesis). 

It has recently been claimed p, ^ that this simple picture may be violated 
if inflation ends with a period of preheating, a violent decay of the inflaton 
particles into another field (or even into quanta of the inflaton field itself). Such 
a phenomenon would completely undermine the usual inflationary picture, and 
indeed the original claim was that large-scale perturbations would be amplified 
into the non-linear regime, placing them in conflict with observations such as 
measurements of microwave background anisotropics. Given the observational 
successes of the standard picture, these claims demand attention. 

In Section ^.12.2| and Ref. | [I05| ], we discuss the general criteria under which 
large-scale curvature perturbations can vary. As has been known for some 
time, this is possible provided there exist large-scale non-adiabatic pressure 



perturbations, as can happen for example in multi-field inflation models |95 



P8| , pOl , pO| , [75| . Under those circumstances a significant effect is possible during 
preheating, though there is nothing special about the preheating era in this 
respect and this effect always needs to be considered in any multi-component 
inflation model. 

In this section we perform an analysis of the simplest preheating model, as 
discussed in Ref. 0. We identify two possible sources of variation of the cur- 
vature perturbation. One comes from large-scale isocurvature perturbations in 
the preheating field into which the inflaton decays; we concur with the recent 
analyses of Jedamzik and Sigl |50] and Ivanov [49| that this effect is negligible 
due to the rapid decay of the background value of the preheating field during 
inflation. However, we also show that in fact a different mechanism gives the 
dominant contribution, which is second-order in the field perturbations com- 
ing from short-wavelength fluctuations in the fields. Nevertheless, we show 
too that this effect is completely negligible, and hence that preheating in this 
model has no significant effect on large-scale curvature perturbations. 
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4.1 Perturbation evolution 



We describe the perturbations via the curvature perturbation on uniform- 
density hypersurfaces, denoted ( In hnear theory the evolution of ( is well 
known, and arises from the non-adiabatic part of the pressure perturbations, 
Eq. ( gl67D . 



On large scales anisotropic stress can be ignored when the matter content 
is entirely in the form of scalar fields, and in its absence the non-adiabatic 
pressure perturbation determines the variation of (, and rewriting Eq. ( p.205| ) 
we get 11, ^ 

7T7 = — OPnad, 4.1 

dN p + P 

where = Ina measures the integrated expansion. The uniform-density hy- 
persurfaces become ill-defined [cf. the definition of (, Eq. (|2.80|) ] if the density 
is not a strictly decreasing function along worldlines between hypersurfaces of 
uniform density, and one might worry that this undermines the above analy- 
sis. However we can equally well derive this evolution equation in terms of the 
density perturbation on spatially-flat hypersurfaces, 6p^ = —{dp/dN)(, which 
remains well-defined. Spatially- flat hypersurfaces are automatically separated 
by a uniform integrated expansion on large scales, so the perturbed continuity 
equation in this gauge takes the particularly simple form 

dSp^ 



dN 



-3(5p^ + 6p^) . (4.2) 



From this one finds that oc dp/dN for adiabatic perturbations and hence 
again we recover a constant value for (. However it is clearly possible for 
entropy perturbations to cause a change in ( on arbitrarily large scales when 
the non-adiabatic pressure perturbation is non-negligible. 

4.2 Preheating 

During inflation, the reheat field into which the inflaton field decays possesses 
quantum fluctuations on large scales just like the inflaton field itself. As these 
perturbations are uncorrelated with those in the inflaton field, the adiabatic 
condition will not be satisfied, and hence there is a possibility that ( might 
vary on large scales. Only direct calculation can demonstrate whether the 
effect might be significant, and we now compute this effect in the simplest 
preheating model, as analyzed in Ref. |]^, |^ . This is a chaotic inflation model 
with scalar field potential 

\/(<^,x) = ^mV + ^^?W, (4.3) 



^^This is the notation of Bardeen, Steinhardt and Turner [Q. General issues of pertur 



bation description and evolution are discussed Section 2.12 and Ref. [105|. The curvature 
perturbation of comoving spatial hypersurfaces, TZ ||6C|, [75[, is the same as — C well outside 
the horizon for adiabatic perturbations, since the two coincide in the large-scale limit, see 



Section 2.12.1 
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where is the inflaton and x the reheat field, and m and g are the infiaton 
mass and the couphng of the reheat field to the infiaton, respectively. Slow-roll 
infiation proceeds with ip > mpi and gx The effective mass of the x fi^ld 

is gip and thus will be much larger than the Hubble rate, H ~ \J Att / mip / m^\, 
for g ^ m/mpi ~ 10~^. Throughout this section, we use the symbol '~' to 
indicate equality within the slow-roll approximation. 

This model gives efficient preheating, since the effective mass of x oscillates 
about zero with large amplitude. In most other models of inflation, preheating 
is less efficient or absent, because the mass oscillates about a nonzero value 
and/or has a small amplitude. 

Any variation of C, during preheating will be driven by the (non-adiabatic 
part of) the x field perturbation. Our calculation takes place in three steps. 
The first is to compute the perturbations in the x field at the end of infia- 
tion. The second is to compute how these perturbations are amplified during 
the preheating epoch by the strong resonance. Finally, the main part of the 
calculation is to compute the change in driven by these x perturbations. 



4.2.1 The initial quantum fluctuation of the x-field 

Perturbations in the x field obey the wave equation 



6x + SH6x + ( ^ + 9'^' ]SX = 0. (4.4) 



The slow-roll conditions ensure that the x field remains in the adiabatic vac- 
uum state for a massive field 



3 



-iujt 



6xk ^ 1=^ , (4.5) 

where tu^ = k'^/a'^+g'^if'^. This is an approximate solution to Eq. ( |4.4| ) provided 

where = gip is the effective mass of the x field. 

Hence the power spectrum (defined in Eq. (|3.55|) ) for long- wavelength fiuc- 
tuations {k -C m^) in the x field simply reduces to the result for massive field 
in fiat space 

Physically, this says that at all times the expansion of the Universe has a 
negligible effect on the modes as compared to the mass. In particular, at the 
end of inflation we can write 
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The power spectrum has a spectral index n^^ = 3. This is the extreme hmit 
of the mechanism used to give a blue tilt in isocurvature inflation scenarios 



|87| 



68. 



4.2.2 Parametric resonance 

After inflation, the inflaton field if oscillates. Strong parametric resonance 
may now occur, amplifying the initial quantum fluctuation in x fo become a 
perturbation of the classical field x- The condition for this is 

g^^»l, (4.9) 
4m^ 

where <^ini is the initial amplitude of the <^-field oscillations. 

We model the effect of preheating on the amplitude of the x fi^ld following 
Ref. 13 as 

nx = nxlend (4.10) 

and the Floquet index is taken as 

/ifc~^ln(l + 2e--''') , (4.11) 

with 

. f * y . ^ f * y . (4.12) 

For strong coupling (g ^ 1), we have -C 1 for all modes outside the Hubble 
scale after inflation ends {k < fcend)- Therefore /ifc ~ ln3/27r 0.17 is only 
very weakly dependent on the wavenumber k. Combining Eqs. ( [4. 8] ) and (4.1(J) 
gives 



4.2.3 Change in the curvature perturbation on large scales 

In order to quantify the effect that parametric growth of the x fi^ld fluctuations 
during preheating might have upon the standard predictions for the spectrum 
of density perturbations after inflation, we need to estimate the change in the 
curvature perturbation C, on super-horizon scales due to entropy perturbations 
on large-scales. 

The density and pressure perturbations due to first-order perturbations in 
the inflaton field on large scales (i.e. neglecting spatial gradient terms) are of 
order g^Lp^x^X- Not only are the field perturbations 5x strongly suppressed 
on large scales at the end of inflation [as shown in our Eq. ( |4.8| )] but so is 
the background field x- We can place an upper bound on the size of the 
background field by noting that in order to have slow-roll chaotic inflation 
(dominated by the m^(y9^/2 potential) when any given mode k which we are 
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interested in crossed outside the horizon, we require x ^ i^/9- The large 
effective mass causes this baclcground field to decay, just like the super-horizon 
perturbations, and at the end of inflation we require x ^ ^/ gik/kcndY^"^ when 
considering preheating in single-field chaotic inflation. Combining this with 
Eq. ( ^.8| ) we find that the spectrum of density or pressure perturbations due to 
linear perturbations in the x fi^ld has an enormous suppression for k <^ fcend^ 

/ \ 3 / \ 2 

477 / m \ ( mpiHc 




( ^ \ ( "^Pl-^cnd \ / k \ 



^x.xU«t^h£- (4.14) 



Effectively the density and pressure perturbations have no term linear in 6x, 
because that term is multiplied by the background field value which is vanish- 
ingly small. 

By contrast the second-order pressure perturbation is of order g'^ip^dx^ 
where the power spectrum of 5x^ is given by |73[ 



/"=-^nx(|k'|)nx(|k-k'|) . , 



27rJo |k|''|k-kf ^ ^ 

We impose the upper limit /ccut to eliminate the ultraviolet divergence 

associated with the vacuum state. Substituting in for Vsx from Eq. ( [4.8|) , we 
can write 

^'''i- = if(5^^) (%^) it^) ih) ■ 

Noting that Hend ~ 'm and /ccut ~ ^max ~ Q^^^kend, it is evident that the 
second-order effect will dominate over the linear term for k < g^^'^q^^^kcnd- 

The leading-order contributions to the pressure and density perturbations 
on large scales are thus 

1 1-2 

6p = m'^Lp5ip + Lp5Lp - Lp'^(l)+ -g^ip^Sx^ + -Sx , (4.17) 

1 1-2 

Sp = -m^ip Sip + p Sip- p^cp - -g'^p^Sx^ + -Sx ■ (4.18) 

We stress that we will still only consider first-order perturbations in the metric 
and total density and pressure, but these include terms to second-order in 6x- 
From Eqs. ( |2.167| ), (|07| ) and ( P^ we obtain 

5pnad = ^ , (4.19) 

where the long-wavelength solutions for vacuum fluctuations in the p field 
obey the adiabatic condition 



^ = . (4.20) 

p p 
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Inserted into Eq. ( [4.1|) , this gives the rate of change of C- 

Note that the non-adiabatic pressure will diverge periodically when = 
as the comoving or uniform density hypersurfaces become ill-defined. Such 
a phenomenon was noted in the single-field context by Finelli and Branden- 
berger |2^, who evaded it by instead using Mukhanov's variable u = aS<f^ 



which renders well-behaved equations. Linear perturbation theory remains 
valid as there are choices of hypersurface, such as the spatially-flat hypersur- 
faces, on which the total pressure perturbation remains finite and small. In 
particular, we can calculate the change in the density perturbation due to the 
non-adiabatic part of the pressure perturbation on spatially-flat hypersurfaces 
from Eq. ( [4.2|) , which yields 

Apnad = -^j SPnadHdt. (4.21) 

Even though 5pnad contains poles whenever = 0, the integrated effect re- 
mains finite whenever the upper and lower limits of the integral are at 7^ 0. 
From this density perturbation calculated in the spatially-fiat gauge one can 
reconstruct the change in the curvature perturbation on uniform density hy- 
persurfaces 

= H^P^ . (4.22) 
Substituting in our expression for 5pnad we obtain 

AC = -^/ {l + ^-^^g'v'\5x'\Hdt, (4.23) 

where we have averaged over short timescale oscillations of the x-field fluctua- 



tions to write 
evolution for t 



. 2 

= \^X^\- To evaluate this we take the usual adiabatic 
le background if field after the end of inflation 

sin(mAt) 

^ = ^ini , (4.24) 

and time-averaged Hubble expansion 

S{mAt + B) 

where G is an integration constant of order unity. The amplitude of the x- 
field fluctuations also decays proportional to 1/ At over a half-oscillation from 
mAt = nvr to mAt = {n + l)7r, with the stochastic growth in particle number 
occurring only when ip = 0. Thus evaluating A( over a half-oscillation At = 
vr/m we can write 

2g^\5x^\xt r"+' f I s\ s\ , 
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where x = mAt, s{x) = sinx/x, Xn = nir and a dash indicates differentiation 
with respect to x. The integral is dominated by the second term in the bracket 
which has a pole of order 3 when s' = 0. Although s/s' diverges, it yields a 
finite contribution to the integral which can be evaluated numerically. For 
Xn ^ ^ the integral is very well approximated by independent of the 

integration constant G. 

This expression gives us the rate of change of the curvature perturbation ( 
due to the pressure of the field fluctuations 6x^ over each half-oscillation of the 
infiaton field ip. Approximating the sum over several oscillations as a smooth 
integral and using Eq. ([4.10|) for the growth of the x-field fluctuations during 
preheating (neglecting the weak /c-dependence of the Floquet index, /x^, on 
super-horizon scales) we obtain 

end _2/tmAt /a o'7^ 

(^na.d — -- 5 '2 ■ 

/TT/i m"^ 

The statistics of these second-order fluctuations are non-Gaussian, being 
a x^-distribution. Both the mean and the variance of Cnad are non-vanishing. 
The mean value will not contribute to density fluctuations, but rather indi- 
cates that the background we are expanding around is unstable as energy is 
systematically drained from the infiaton field. We are interested in the variance 
of the curvature perturbation, and in particular the change of the curvature 
perturbation power spectrum on super-horizon scales which is negligible if the 
power spectrum of Cnad on those scales is much less than that of ( generated 
during inflation, the latter being required to be of order 10~^° to explain the 
COBE observations. 

To evaluate the power spectrum for Cnad we must evaluate the power spec- 
trum of Sx^ which is given by substituting Vsx, from Eq. (|4.13| ), into Eq. (|4.15| ). 
This gives 



where 



I{K,mAt) = - c/k' / rf^e2(^-'+'^--')"^*/€'2sin^, (4.29) 
2 Jo Jo 



K = fc/fcmax as defined in Eq. (|4.12|) , and 6* is the angle between k and k'. Note 
that at the end of inflation we have /(k,0) = k^^^ ~ 1, and Vs^i oc k^. This 
yields 



One might have thought that the dominant contribution to Cnad on large 
scales would come from Sx fluctuations on those scales, and that is indeed the 
presumption of the calculation of Bassett et al. ^. However, in fact the 
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-20 -15 -10 -5 

Figure 5: The power spectrum of the non-adiabatic curvature perturbation 
V(^^_^^, shown at four different times: from bottom to top niAt = 0, 50, 100 
and 150. The parameters used were g = 10^^, m = lO^^mpi and kcut = k^ax- 



integral is initially dominated by k' ~ k^nt, namely the shortest scales. The 
reason for this is the steep slope of Vsx] were it much shallower (spectral index 
less than 3/2), then the dominant contribution would come from large scales. 

To study the scale dependence of J(k, mAt) and hence "P^^^^ at later times, 
we can expand fi^-K' for kk' <^ 1 as 

2k' cos ^ 9x N 

= + ^-^^ + OiK^) . (4.31) 



We can then write the integral in Eq. ( [4.29| ) as 

mAt) = /o(mAt) + 0{k^) , (4.32) 
where first-order terms, 0{k,), vanish by symmetry and 

/o(mAt) = - e^^'=""^*/€'2rffi:' . (4.33) 

2 Jo 

Thus the scale dependence of Vq^^ remains k^ on large-scales for which k ^ 1. 

At late times these integrals become dominated by the modes with k'"^ <^ 
{mAt)^^ which are preferentially amplified during preheating. These are 
longer wavelength than kcnt, but still very short compared to the scales which 
give rise to large scale structure in the present Universe. From Eq. ( 4.11 ) we 
have /ifj/ ~ /io — '^'^/3, for k'^ -C 1, where /ig = (ln3)/27r, which gives the 
asymptotic behaviour at late times 

Jo ~ 0.86(mAt)-=^/2g4MomAi _ ^4 34^) 
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Thus although the rate of growth of V(^^^ becomes determined by the expo- 
nential growth of the long-wavelength modes, the scale dependence on super- 
horizon scales remains proportional to for k < (mAt)"^/^. This ensures 
that there can be no significant change in the curvature perturbation, (, on 
very large scales before back-reaction on smaller scales becomes important and 
this phase of preheating ends when mAt ~ 100 . 



Numerical evaluation of Eq. ( ^.301 ) confirms our analytical results, as shown 
in Fig. ^ For k <^ femax, the spectral index remains k^ during preheating. Ob- 
servable scales have log^g k/kmux — —20. 



4.3 Conclusions 

Our result shows that because of the k^ spectrum of 6x, which leads to a 
similarly steep spectrum for ^nad, there is a negligible effect on the large-scale 
perturbations before the resonance ceases. The fluctuations in x grow largest 
on small scales, so that backreaction can turn off the resonant amplification 
on these scales before the curvature perturbation can become non-linear on 
very large scales. 

This contrasts with the original findings of Ref. [^, in which Bassett et 
al. report the results of their investigation of the model given by Eq. ( |4.3| ). 
They solved the field equations for particular modes to linear order numerically, 
assuming that the initial x spectrum is flat. For strong parametric resonance 
this would lead to an exponential amplification of the curvature perturbation 
( on super-horizon scales into the non-linear regime before backreaction can 
shut down the resonance. This would undermine the standard cosmological 
model, in which the curvature perturbation is constant on large scales. 

Our result was derived analytically, with exemption of the integral in 
Eq. ( |4.26| ), which we had to evaluate numerically. The use of an analytic 
approximation, as described in Section ^]2|, allowed us to calculate the power 
spectrum of the x fluctuations and hence enables us to calculate their effect 
over the whole range of scales. The suppression of the large-scale perturba- 



tions in 5x, discussed in Refs. I^O], ^], means that large-scale perturbations 



in 6x are completely unimportant. However, it turns out that they don't give 
the largest effect, which comes from the short-scale modes which dominate the 
integral for (nad- Nevertheless, even they give a negligible effect, again with a 
k^ spectrum. Indeed, that result with hindsight can be seen as inevitable; it 



has long been known ||108i pq , |T6| that local processes conserving energy and 
momentum cannot generate a tail shallower than k^ (with our spectral index 
convention) to large scales, which is the Fourier equivalent of realizing that in 
real space there is an upper limit to how far energy can be transported. Any 
mechanism that relies on short-scale phenomena, rather than acting on pre- 
existing large-scale perturbations, is doomed to be negligible on large scales. 
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5 Concluding remarks and discussion 



5.1 Summary 

In Section ^ we presented the framework of cosmological perturbation theory 
for single and multiple fluids. This allowed us to introduce the notion of 
gauge-invariant perturbations and to provide their governing equations, which 
we then used throughout the rest of this work. 

In Sections f^.l3| and we have identified the general condition under 



which the super-horizon curvature perturbation on spatial hypersurfaces can 
vary as being due to differences in the integrated expansion along different 
worldlines between hypersurfaces. As long as linear perturbation theory is 
valid, then, when spatial gradients of the perturbations are negligible, such a 
situation can be described using the separate universes picture, where regions 
are evolved according to the homogeneous equations of motion. 

Equivalently, we can study the model under consideration using the per- 
turbed energy conservation equation. As discussed in detail in Section |2.12 



and in Ref. ||1U5 |, large-scale curvature perturbations can vary provided there 
is a significant non-adiabatic pressure perturbation. This is always possible in 
principle if there is more than one field or fluid, and since for example pre- 
heating usually involves at least one additional field into which the infiaton 
resonantly decays, such variation is in principle possible during preheating. 
In particular, the curvature perturbation on uniform-density hypersurfaces, 
C, can vary only in the presence of a significant non-adiabatic pressure per- 
turbation. The result follows directly from the local conservation of energy- 
momentum and is independent of the gravitational field equations. Thus ( is 
conserved on sufficiently large scales in any metric theory of gravity, including 
scalar-tensor theories of gravity or induced four-dimensional gravity in the 



brane-world scenario |]7S|, |105| . 



Multi-component infiaton models are an example where non-adiabatic per- 
turbations may cause the curvature perturbation to evolve on super-horizon 
scales. 



In Section |^ we have shown that the recently proposed model of assisted 
inflation, driven by many scalar fields with steep exponential potentials, can be 
better understood by performing a global rotation in field space, which allowed 
us to re-write the potential as a product of a single exponential potential for 
a weighted mean field, (f, and a potential Vn for the orthogonal degrees of 
freedom, cxj, which has a global minimum when cxj = 0. 

We have also been able to give exact solutions for inhomogeneous linear 
perturbations about the attractor trajectory in terms of our rotated fields. 
Perturbations in the weighted mean field ip corresponds to the perturbations 
in the density on the uniform curvature hypersurfaces, or equivalently, pertur- 
bations in the curvature of constant density hypersurfaces and hence constitute 
perturbations along the attractor trajectory. They therefore describe adiabatic 
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curvature perturbations. 

First-order perturbations in the fields orthogonal to the weighted mean field 
are isocurvature perturbations during infiation. Vacuum fiuctuations lead to 
a positively tilted spectrum. We have shown that the non-adiabatic pertur- 
bations decay relative to the adiabatic perturbations and hence we recover 
the single field limit at late times. In particular we find that the curvature 
perturbation ( becomes constant on super-horizon scales during infiation. 

Recently field rotations have been used in a more general context than 
assisted infiation [^]. It has been shown that for multiple scalar fields a lo- 
cal rotation in field space can always be used to separate out adiabatic and 
entropic modes. It was shown that the non-adiabatic part of the pressure 
perturbation on large scales is proportional to the curvature of the trajectory 
in field space, and that ( remains constant on large scales for straight line 
trajectories such as the assisted infiation attr actor with (t^ = 0. 

In Section ^ we have focussed on the simplest preheating model, as dis- 
cussed in Ref. 0. We have identified the non-adiabatic pressure, and shown 
that the dominant effect comes from second-order perturbations in the pre- 
heating field. Further, the effect is dominated by perturbations on short scales, 
rather than from the resonant amplification of non-adiabatic perturbations on 
the large astrophysical scales. Nevertheless, we have shown that the contri- 
bution has a spectrum to large scales, rendering it totally negligible on 
scales relevant for structure formation in our present Universe by the time 
backreaction ends the resonance. Amongst models of infiation involving a 
single-component infiaton field, this model gives the most preheating, and so 
this negative conclusion will apply to all such models. 

In Ref. [0 Bassett et al. have suggested large effects might be possible in 
more complicated models. They consider two types of model. In one kind, 
inflation takes place along a steep-sided valley, which lies mainly along the 
direction of a fleld cp but with a small component along another direction x- 
In this case, one can simply deflne the inflaton to be the fleld evolving along 
the valley floor, and the second heavy fleld lies orthogonal to it. Taking that 
view, there is no reason to expect the preheating of the heavy fleld to give rise 
to a bigger effect than in the simpler model considered in this paper. 

In the second kind of model, the reheat fleld is light during inflation, and 
this corresponds to a two-component inflaton fleld. As has long been known, 
there can indeed be a large variation of ( in this case, which can continue 
until a thermalized radiation-dominated universe has been established. In- 
deed, in models where one of the flelds survives to the present Universe (for 
example becoming the cold dark matter), variation in ( can continue right to 
the present. This variation is due to the presence on large scales of classical 
perturbations in both flelds (properly thought of as a multi-component infla- 
ton fleld) generated during inflation, and the effect of these must always be 
considered in a mult i- component inflation model, with or without preheating. 
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Recently, using the results presented in Sections and H, we were able to 
present the first calculations predicting the over-production of primordial black 



holes (PBHs) due to preheating [33|. The scales which we are interested in 
pass outside the Hubble radius towards the end of inflation. The fluctuations 
in the preheat field x ^^re amplified during preheating and large PBHs will be 
formed, when these scales re-enter the Hubble radius after reheating during 
the radiation dominated era, if the fluctuations are sufficiently large. We found 
that for a wide range of parameter values in the massive infiaton model given 
by Eq. (|4.3|) PBHs are over-produced before back-reaction can shut down the 



resonant amplification. This poses a serious problem for the standard theory 



of preheating as described in |57 



5.2 Extensions 

A natural extension of the cosmological perturbation theory of 4-dimensional 
spacetimes, perturbed around a FRW background, is to higher dimensional 
theories. Non-Einstein gravity (in our four spacetime dimensions) may for 
example emerge from theories involving a large extra dimension [^, 



p8| , p9| . Recently there has been a lot of interest in this so called brane- world 
scenario In this model matter is confined to the world volume of a 4- 
dimensional three-brane, whereas gravity lives as well in the higher dimensional 
bulk spacetime. The simplest case of such a higher dimensional spacetime is 
a 5-dimensional spacetime, in which 4-dimensional slices correspond to "our" 
standard FRW universe. 

In order to study inhomogeneous perturbations we will pick a specific form 
for the unperturbed 5-d spacetime metric that accommodates spatially fiat 
FRW cosmological solutions on the brane, 

dsl = —n'^{t, y)drf' + a^(t, y)5ijdx^dx^ + dy^ , (5.1) 

where a = a{t, y) and n = n{t, y) are scale factors and y is the "new" fifth 
dimension. The line element given above includes anti-de Sitter spacetime as 
a special case and has been extensively studied in the literature 0. 

We consider arbitrary linear perturbations about the background metric. 
In keeping with the standard approach of cosmological perturbation theory 0, 
as presented in Section ^ of this thesis, we will introduce scalar, vector and 
tensor perturbations defined in terms of their properties on the 3-spaces at 
fixed t and y coordinates. 

We can write the most general metric perturbation to first-order as 

/ -n2(l + 20) na{B\i - Si) n(j)y 

9ab = na{By - Sj) [{1 - 2tlj)-fij + 2E\ij + 2F(i\j) + hij\ a{By\i - Syi) 

\ nCpy a{By\i - Syi) 1 + 2(t)yy 

(5.2) 

where 0, B, ip, E, (j)y and (j)yy are scalar s, 5, F, and Sy are (divergence- 
free) 3-vectors, and hij is a (transverse and traceless) 3-tensor. The reason 



63 



for splitting the metric perturbation into these three types is that they are 
decoupled in the linear perturbation equations, as in the 4-dimensional case 
(for a separable metric Ansatz). 

In the perturbed spacetime there is a gauge-freedom in the definition of the 
scalar and vector perturbations, as in the standard case described in Section 0. 
Under a first-order coordinate transformation, ^ + which we will 
write here as 

t t + 6t 
x' x' + 6x^' + 6x' (5.3) 

y y + 5y 

where St, Sx and Sy are scalars and 6x^ is a (divergence-free) 3-vector, the 
perturbations transform as 

fl 111 Tl d ' 

(/) (f)-6t 6t Sy , B ^ B + -St Sx , 

n n an 

Si — >■ Si -\ — Sxi , By ^ By ~ aSx' Sy , 

n ay ^ 

Qj q! 

Syi Syi + aSxi , ^ %Ij + -St H Sy , (5.4) 

(X ct 

E E-Sx, F,-^ Fi-Sxi, 

(j)y (f)y + nSt' Sy , (f)yy (f)yy - Sy' , 

where dot and prime denote differentiation with respect to coordinate time t 
and bulk coordinate y, respectively. Comparing these equations with the trans- 
formation properties of scalar and vector perturbations in the 4-dimensional 
case, Eq. ( p.21|) -( p.26|) , we see that there is now even more scope for "gauge- 



ambiguities". But we can easily circumvent any such difficulties by, again, 
simply fixing the gauge and constructing gauge-invariant quantities. For ex- 
ample, we can construct two gauge-invariant vector perturbations from the 
vector perturbations given above, 

~S, = S, + -Fi, (5.5) 
n 

Syi = Syi + aFl . (5.6) 

Although we have "only" introduced one more dimension to the problem 
the number of degrees of freedom increases quite dramatically. In the 4- 
dimensional case we had four scalar perturbations, 0, B, ip and E, and two vec- 
tor perturbations, Fj and Si. This allowed us to construct two gauge- invariant 
scalar perturbations and one vector gauge-invariant perturbation (see Section 



2.5) 



In the 5-dimensional case we have seven scalar perturbations, 0, B, ip, E, 
By, (py and (pyy, sjid three vector perturbations, Fi, Si and Syi. This allows us 
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to construct four gauge-invariant scalar perturbations and two vector gauge- 
invariant perturbations. The theory of cosmological perturbations in more 
than four dimensions will therefore be much richer than in four dimensions. 

First steps towards solving these problems have already been undertaken 
p5| , p^ , p9| , p9| . For example, our proof of the constancy of the curvature 
perturbation ( on large scales, Eq. (|2.204| ), which follows from the local con- 
servation of energy-momentum independent of the form of the gravitational 



field equations, validates a recent discussion |]78i of chaotic inflation in the 



brane-world scenario, which relied on that equation. 
5.3 Outlook 

The theory of cosmological perturbations has made dramatic progress since 
the seminal works of Lifshitz and Bardeen [0 . Nevertheless, the theory is 



far from complete with many open questions remaining. 

What should be done next? Although the governing equations for multi- 



component systems with energy exchange have been developed and applied 
to simple systems of fluids quite some time ago, the system of equations 
used is not entirely satisfactory. Using the formalism outlined in this work in 
Section ^ and using insight into mult i- component systems gained in Section 
and it should be possible to cast the governing equations for multi- 
component systems with energy exchange into a new gauge-invariant form. 
Splitting the system again into adiabatic and entropic perturbations it should 
be possible to show that the adiabatic perturbations are sourced by the en- 
tropic perturbations, whereas the entropic perturbations can't be sourced by 
the adiabatic ones on large scales. Interesting systems would, for example, 
include quintessence or a form of coUisional dark matter as well as "standard" 
forms of matter. 

One obvious extension of the cosmological perturbation formalism is to 



higher dimensional theories, as outlined in the previous section |5.2| . Even if it 
should turn out that we do not live on a 4-dimensional 3-brane embedded in 
a higher dimensional space, we can learn a lot about standard 4-dimensional 
cosmological perturbations and improve our technical expertise from investi- 
gating these theories. 

Recently we were able to show in Ref. ||3^ that the preheating model stud- 



ied in Section |^ leads to an overproduction of primordial black holes (PBH) 
during the radiation dominated era before backreaction can shut down the 
amplification of the field fluctuations. This constitutes the first consistent 
investigation of this kind. But is this a generic feature of preheating, or is 
there something wrong with the backreaction mechanism? It would therefore 
be interesting to study other models of preheating with respect to their PBH 
generation behaviour. 
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A Appendix 
A.l Notation 



The sign convention is (+++) in the classification of Ref. ||82[ . 

The number of spatial dimensions is throughout the main part of the thesis 

n = 3, but we do not fix n in the appendix. 

Tensor indices: 

Greek indices, such as a, /5, . . . , /i, i/, . . ., run from to n, that is over all 
dimensions. Latin indices, such as a,b, . . . ,i, j, . . ., run from 1 to n, that 
is only over spatial dimensions. 
The connection coefficient is defined as 

= ^9"^' {9c.P,^^ + ^"M,/3 - 9Pt.,a) ■ (A.l) 
The Riemann tensor is defined as 

pa -pa po I -pa. pA -pa -p\ ( \ 0\ 

The Ricci tensor is a contraction of the Riemann tensor and given by 

and the Ricci scalar is given by contracting the Ricci tensor 

R = R''^. (A.4) 

The Einstein tensor is defined as 

Gij,y = R^u — ij^9^iuR ■ (A-5) 

The covariant derivatives are denoted by 

= Covariant differentiation with respect to g^v , 
|j = Vj Covariant differentiation with respect to 'jij . 

Throughout this work we use the units c = h = 1. 
A. 1.1 Time variables 

Throughout this work we use as time variable the conformal time 1], if not 
stated otherwise. Conformal time is related to coordinate time t by 

dt = adt] . (A. 6) 

The time derivatives are denoted by 

= — Time derivative with respect to coordinate time , 
dt 

' = — Time derivative with respect to conformal time . 
dr] 



A.l 



Another useful time variable is proper time r. The definition follows from the 
line element ds^ in Eq. (|2.11|) as 



dr = a^/l + 2(f) dr] , (A. 7) 

and can be rewritten, to first-order in 0, as 

dr = a {1 + (p) dr] . (A.8) 

A. 1.2 Metric perturbation variables 

Before relating our notation to other works, we give the linearly perturbed 
metric tensor around a FRW background already presented in Section |2.1| . 
The covariant metric tensor is given by 



and the contravariant metric tensor is 
9"' = a-\v) 



-(1-20) B^'-S' 
fi| ^ - (1 + 2ij)Y^ - 2E| - F] ^ - F\ ' - W 



(A.IO) 

The notation of Kodama and Sasaki [55| is related to our notation as follows 
The scalar perturbations are 

E^^, ^P^-(^Hl + Iht^Y. (A.ll) 

Another useful relation is ■0 — ^^"^E = —HlY. The alternative notation 



7?. = — ^ is also used. The vector perturbations are 



Y 



S,^B^'^y}'\ = (A.12) 



The tensor perturbations are related by 



The scalar harmonic functions are defined as 

V^F = -k^Y, Yi = -^Y. 

k 

1 



K,=2HPY^f . (A.13) 



1 1 



Y,, ^ + (A- 14) 



^^The notation of Bardeen [g| is similar to that in but the harmonic functions are 
denoted Q(o). 



A.2 



where = V*Vi. The vector valued harmonic functions are defined by 

(v= + F)fW = 0, f'"* = o, 

Note as well that FAj-* = ^ {k'^ — {n — 1)k) Y-^^\ The tensor valued harmonic 
functions are defined by 

(V2 + A;2)Ff = 0, Y^'^\ = 0, 

Y^^% = 0. (A. 16) 



A. 1.3 Matter variables 



The scalar matter variables in the notation used by us compare to the variables 
used by Kodama and Sasaki in |E5| as 



k ' 

Sp^pttlY, Tr = ^7rTY, (A. 17) 

and n^j = pTixY^y The vector matter variables are related by 

t;, = t;«F,, 7^i,\,-)=4^Yi^\ (A.18) 
and the tensor matter variables by 

(*™)7rij- = . (A. 19) 

A. 2 The connection coefficients 

The connection coefficients in the FRW background are 

r° — h r° — n 

r}o =h6]. r5, = (")r},, (A.20) 

where '•"■'r*^ is the connection of the n-sphere. The perturbed connection 
coefficients for the scalar perturbations corresponding to the metric given in 
Eq. ( ^lOD are 



^r^o =0', r°, = 0, + /is,, 



5V% =-(B-E'- 2hE\^^ + {-2h (0 + ^) + ^') li, , 
5Vl, = [hB + B' + ct>)\\ 

A.3 



The perturbed connection coefficients for the vector perturbations are 



2 



For the tensor perturbations we find 



00 



pO 



Oi 



00 



0, 



Oj — 2 0' 



5r 



jk 



-h[j + h hij , 



(A.22) 



(A.23) 



Note that h is the Hubble parameter in conformal time and not related to the 
tensor perturbation hy. 



A. 3 The Einstein tensor components 

In the background the Einstein tensor takes the form 

n{n — 1) 



2a2 



{h' + k) , 



G'. = ri-l ^a" ^ n-A 



7 2 n — 2 



For the scalar perturbations we find 



Note that Gl - nG\ = {h^ - h' + k). 



5G^ 



5G^j 



n — 


1 


a2 




n — 


1 


a2 




n — 


1 


a? 



(A.24) 



[h<l> + ^' + K{B-E%,, 



2— + {n- 4)/i2 ^(f) + h(l)' + < + hip+{n- 2)hip 



— (n — 2)ft;V' 



a 
A.4 



(A.25) 



where D = (p — {n — 2)ip — ha — a' — {n — 2)ha. The trace of the spatial part 
of the perturbed Einstein tensor is obtained by contracting Eq. ( |A.25| ) to get 



n{n — 1) 



2— + in- 4:)h^ [(p + h(f)' + ip" + h-^ + (n - 2)hip 
a 



— {n — 2) Kip 



+ 



(n- 1) .r. 



(A.26) 



For the vector perturbations we find 



6G\ 
5G\ 



{n - 1)k + 
2^2 



(A.27) 



Note that the last expression in Eq. ( A.27 ) above can be rewritten in terms of 
the vector shear Tij, given in Eq. (|2.42| ), as 



6Gij 



[n 



2)hn,] . 



(A.28) 



For the tensor perturbations the Einstein tensor components are 



6G' 



0, 



= 



A.5 



A.4 List of Symbols 



a 


Scale factor 




Acceleration 


Cs 


Adiabatic sound speed 




Adiabatic sound speed of Jth component 


f(J) 


Momentum transfer perturbation of Jth component 


g 


Coupling of preheat field x to the inflaton 


g^v 


Metric tensor 


gij.v 


Conformally rescaled metric tensor 


h 


a' /a 


hij 


Tensor metric perturbation 


k 


Comoving wavenumber 


m 


Mass of the inflaton field 


mpi 


Planck mass 




Mass of preheat field x 


n 


Scale factor 


ris 


Spectral index of curvature perturbations 


Anx 


Scale dependence of perturbation spectrum of a quantity X 


P 


Pressure 


Pj 


Pressure of Jth component 


Pi 


Exponent of ith exponential potential 


P 


"Mean" exponent 


Q 


Dimensionless coupling parameter 


Q(J) 


Reduced energy tranfer parameter of Jth component 


ds 


Infinitesimal line element 


t 


Coordinate time 


At 


Coordinate time elapsed since end of inflation 


u 


Rescaled perturbation in mean scalar field, aSip 


u^' 


4-velocity 


V 


Scalar velocity perturbation 


Vj 


Scalar velocity perturbation of Jth component 


Vj 


Rescaled perturbation in Jth orthogonal field, aSaj 


v' 


Vector velocity perturbation 




Spatial coordinate 


y 


"Fifth" dimension 


z 


Mukhanov potential, a'^(p'/a' 



A.6 



R 
u 


OllliL VCOLUl I oOdidi lliCLlll^ ^kjL LUi UdLlUli J 






C 


ScriliriP' solntion 

v^CAi XX £1, kJ V_y X iX U X V^XX 


E 


Anisotropic stress perturbation (scalar metric perturbation) 




Vector metric perturbation 


G 


Newton's constant 


Ljr 111/ 


Einstein tensor 


H 


Hubble parameter with respect to coordinate time 


^(1) 


rn Ti uf^i Til Ti r*"!"! on "f" ri <^ Ti T'G'l' In TiH r\T rlf^cfirf^f^ 7/ 

iidlixvCl i LiliL. ulwli \JL uliC' iii O L JVlllU, Ul ^J.Ccii KjKj ly 


K 


FiXtrinsir rnrvr^turp 

X i./V U X X XXlJ X LXX V CXj U l^X V> 




J^CXjCIX CXjXXCI XCXjXX 


c 


LaEfranEfian density 

J-JCAj^^X CAjXX Li^X CAjXX V_X^^XXk_JX u y 


N 


Number of e-folds (integrated expansion) 




Unit time-like vector field 


p 


Projection tensor, Qp^y -\- N^N^ 


V(J) 


Energy transfer parameter of the Jth component 


Q{j) 


Energy momentum four vector of the Jth component 


rx 


l-^/'^^irOT' GT\£i/^'l'T'n TYI i^T £1 /"111 Ci TI'I'tI'TT' X 

jruwtJi opccti Liiii ui <x mj-diitity y\. 


p? 

± L 










Intrinsic snatial curvature 

J-XXUX XXXk^XV.^ L_/CAj UXCAjX ^_4.X V CAJ <J 1_4.X \^ 


n 


Curvature perturbation in comoving gauge 


s 


Action 

X ^v^uxvyxx 




Vector metric perturbation 


s ■ 


Vector metric nertiirbation i braneworld 1 

V U V_y X XXXV_j U X X V-f V-'X U LXX F .7 Ca/ UX V-/XX I P^X C^XXVy V V V-/X X\_X 1 


Stt 


Entropy perturbation 


Sij 


Reduced entropy perturbation 


T 


Energy momentum tensor 


V 


Velocity perturbation in comoving total matter gauge 




Potential of scalar field 


% 


"Mean" potential of n scalar fields 
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X 




AM 




C / -r\ 




V 


V^'UlliUl lllfil Lllllc 


s 




lik 


iViclllC IcllbUl Ull O-J--' bpdCc Willi CUllbldlll CUlVdlLllc 1% 


1^ 

nj 


U.1 Vclt 111 tJ Ul Ud^^lv^lUllllU. bJJdCCLllllC 




r lUC^Util lllU.t;-X- 


LO 


UlIIKsilblOIllGbb ordllb-J-'lCKc pd;! dmytcr 




voiticity 


TT 


/\lllbO II OpiC b II ebb VcClOI 




/\IllbU II UpiC b II ebb IcIlbUI 


tensor „i 


iClioVjl Idl dliloVJ tlijpi*^ totlCtoO tCllOiJl 




OCclldl ilclQ 


<^ 


VVclgiliLcQ llicdli iiclQ 


<^dil 


LyllilcIiblUIiltjbb LllldlUIi iitjlU. 




Tni "^1 Cil Q Tn"r\l I'J'nri o r\T mfi Q'J'j^ti j^g/^iI 1 Q'J'ij^ti 
lllltldl dllipil t UU-t? Ul llllldtUll UbCllldtlUll 


rh 


ijdpbt; lUIlCllUIl ^^bCdldl Iilt:lIlC pel I Ul UdllUIl J 


<Py 


Extra sccilcir potenticil in braneworld 


'ryy 


Il/Xtld bCdldl pOtclKldl 111 DIdllcWOilQ 


tp 


Ul Vd I Ul (3 ptil lUI UdllUIl l^bCdldl IlltilllC ptii L Ul UdLIOIlJ 


P 


n .Tl/^TCnT" H^^TIGtI'IT' 

J-JllCl&y iJ-dloltV 


HJ 


n ,T1 i^T'0"\T H f^Tl CItI'TT' i^T Tl" Tl /^l^TTlT^i^Tl ATI"!" 
J-JliCl^j' U-Cliolty Xjl. U til *^VJllipiJlldlt 


(7 


Olicdi bCdldi 


Cj 


J 111 ilLllI Ul LliU^UIidl tU Illcdli litjlU. 




Ollcdl LcllbUl 


Tij 


VjrdUgc IllVdIIdlll VtJClUI ptil lUI UdllUIl 




r ,"VT\?1 Tl *31 OTl 
J—JjY L/dllOlUll 


S 


Arl~)itrarv srfilrir fiinrtion 


s 


ArV)itrarv sralar fiinrtion 

X XX P^X U X CAiX y l_j v^CAjXCAjX X l_XXX U X V./ XX 


s 


A T r»Tt"Tfi r\T \r(^r*i' f^T ^rsi 1 n rn n r*"!" i/^ti 

iT-lUlLldlj VCL^LUi VdlUCU. ILlliL^LlUli 


A 


DlTTlPTlSlOTllpSS dpTlSltv DPrtlirllRtlOTl 

J— ' xxxxvyxxoxvy xxxvyOkj vxvyxxkjx u y Lyv^x u uix i^cvuxwxx 


^ J 


Dimftnsionlftss dftnsitv Dftrturbation of Jth coniDonftnt 

J — y xxxxv/xxk^xvy xxx\_>k_JkJ v^vyxxkjx u y vyx u lax ca/ ux vyxx v-/x f.j uxx xxx v-/xx\_yxx u 


A T 


r? n /^firi H 1 TYi on ci on 1 agg H ^^tigtI'tt' T^f^T'l'nT nCi'l'n^Ti r^T T"!" n /'T'^tti t^T'^ti at\1~ 
ivtouiiut/U. Liiiiidioiuiiic&o udi&itv pci t U.1 UdtlUll Ul U til *^uiiipuiidit 




H Tl 4" T/^^T^T r T\ OTT n T r~» Q T 1 /~\T1 

H/lllIUpy pel LUI UdLlUIl 


i \X) 


1 -1 Q T > 1 T > 1 Q Tl 1 Tl f^T 1 r~\Tl 

Vjrdlllllld lUliL^LlUli 


J- int 


TTi'l"nTi*3iP PTi'i'TO'nAA TiPT'i'nT'nPi'l'inTt 
xiitiiiioi*^ ciiti ujjy ijci t Lii L/dtiuix 


Trel 


Relative entropy perturbation 


rj 


Intrinsic entropy perturbation of Jth component 


n 


Scalar anisotropic stress tensor 




Bardeen potential (lapse function in longitudinal gauge) 




Bardeen potential (curvature perturbation in longitudinal gauge) 
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